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ABSTRACT 

This  thesis  examines  the  use  of  heat  for  producing  clinical  lesions 
in  tissue  and  presents  the  design  and  analysis  of  a  resistively  heated 
surgical  probe.  The  probe  temperature  is  accurately  maintained  and 
measured  by  using  a  Wheatstone  bridge.  The  probe  was  embedded  in  an 
agar-water  test  medium  and  the  temperature  field  emanating  from  the 
probe  was  measured  using  liquid  crystal  tapes.  The  experimental  results 
compare  within  approximately  10%  of  a  two-dimensional  numerical  solution 
A  one-dimensional  theoretical  model  is  developed  and  results  are 
presented  to  show  the  relative  effects  of  metabolism  and  blood  flow  on 
the  temperature  field. 
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NOTATION 


Symbol 

Description 

Typical  Units 

a 

Tissue  thermal  diffusivity 

2 
cm  /sec 

a' 

Temperature  coefficient  of  resistance 

ohm/ohm/°C 

c 

Tissue  specific  heat 

cal/gm/°C 

cb 

Blood  specific  heat 

cal/gm/°C 

I 

Current 

amps 

k 

Tissue  thermal  conductivity 

-  cal/cm/°C/sec 

mb 

Volumetric  blood  flowrate 

gm/cm  /sec 

Q 

Heating  rate 

watts 

r 

Radial  location  in  tissue 

cm 

ro 

Probe  radius 

cm 

R 

Resistance 

ohms 

m 

Matabolic  heat  generation  rate 

3 

cal/cm  /sec 

t 

Time 

sec 

T 

Tissue  temperature 

°C 

Tb 

Arterial  temperature 

°C 

T 
c 

Critical  temperature 

°C 

To 

Initial  tissue  temperature 

°C 

TP 

Probe  temperature 

°C 

V 

Voltage 

volts 

p 

Tissue  density 

gm/cm 

p* 

Resistivity 

ohm-cm 

Non-dimensional  Groups 

2 
b  b  o 
k 

-4              -: 
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I.  INTRODUCTION 

A  variety  of  surgical  techniques  are  presently  used  to  produce 
clinical  lesions.  In  addition  to  the  familiar  mechanical  or  "scalpel" 
method,  surgical  procedures  employing  chemical,  ultrasonic,  focused 
x-rays,  radiation,  and  cooling  and  heating  agencies  are  used  to  destroy 
descrete  regions  of  tissue,  normally  nervous  tissue.  The  choice  of 
methods  depends  largely  on  the  nature  of  the  clinical  problem.  Ideally, 
the  surgeon  desires  a  technique  that  possesses  the  following  character- 
istics: 1)  Safety,  2)  Reversibility,  3)  Reproducibility,  4)  Sharp 
delimitation,  5)  Hemostasis,  6)  Flexibility,  7)  Simplicity,  and  8) 
Speed  of  application.  At  times,  it  is  not  possible  to  meet  all  of  the 
above  surgical  ideals.  In  such  cases,  the  choice  of  procedures  is 
based  largely  on  an  optimization  of  the  ideals,  with  safety  being  of 
paramount  importance.  Gengler  [10]  presents  a  summary  of  the  various 
methods,  complete  with  the  advantages  and  disadvantages  of  each. 

In  this  thesis,  the  use  of  heat  for  producing  lesions  is  investigated. 
In  particular,  the  design  of  a  resistively  heated  surgical  electrode, 
henceforth  referred  to  as  a  probe,  is  presented.  Tissue  destruction  is 
accomplished  by  using  the  probe  as  a  heating  element  to  raise  the 
temperature  of  a  predetermined  region  of  tissue  above  a  critical  level, 
normally  taken  to  be  55°C  [2]. 

The  unique  feature  of  this  probe  is  that  its  surface  temperature 
can  be  accurately  measured  and  controlled  by  using  a  Wheatstone  bridge. 
As  a  result,  the  temperature  field  in  the  medium  surrounding  the  probe 
can  be  analytically  predicted  by  solving  the  heat  equation  prior  to 
applying  the  probe.  This  is  quite  significant  when  one  considers  that 
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the  volume  of  tissue  destroyed  depends  solely  on  the  local  tissue 
temperature. 

Temperature  fields  determined  experimentally  with  the  probe 
embedded  in  a  0.30%  agar  -  99.7%  water  test  medium  compare  within  an 
accuracy  of  13%  of  those  predicted  using  the  analytical  model.  A  novel, 
relatively  new,  temperature  sensing  device  was  used  in  the  experimental 
phase  of  the  study.  Liquid  Crystals  [12],  a  material  that  changes  color 
over  a  known  and  well  defined  temperature  range,  were  employed  to 
obtain  a  two-dimensional,  visual  display  of  the  transient  temperature 
field  which  developed  around  the  probe.  By  using  two  separate  orien- 
tations of  the  liquid  crystals,  a  three-dimensional  picture  of  the 
temperature  field  was  inferred.  In  addition  to  the  highly  desirable 
visual  aspect  of  the  liquid  crystal  material,  it  also  produces  a  minimum 
disturbance  in  the  test  medium. 

The  effects  of  blood  flow  and  metabolism  on  the  size  of  the  lesions 
produced  were  also  investigated  using  an   analytical  model. 
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II.  SUMMARY  OF  HEATING  TECHNIQUES 

A.   DIRECT  CURRENT  (dc)  METHOD 

Localized  lesions  have  been  produced  in  tissue  by  applying  a  direct 
current  between  two  electrodes.  An  active  electrode  is  inserted  into  a 
predetermined  location  in  the  tissue,  and  an  inactive  electrode,  or 
ground,  is  placed  on  some  other  part  of  the  body.  The  resistivity  of 
the  tissue  causes  electrical  energy  to  be  converted  into  heat,  thus 
increasing  the  local  tissue  temperature.  When  the  tissue  temperature 
is  raised  above  55°C  [2],  the  tissue  is  destroyed. 

Although  small  and  discrete  lesions  can  be  produced  by  the  dc 
method,  it  has  several  disadvantages.  Direct  current  follows  the  path 
of  least  resistivity  and  tends  to  have  preferred  paths  in  tissue,  thus 
causing  irregular  lesions.  Gas  bubbles,  a  result  of  electrolysis,  form 
around  the  active  electrode  and  block  the  flow  of  current.  In  addition, 
dc  devices  suffer  from  polarization  which  can  produce  transients  that 
stimulate  muscle  and  nervous  tissue  [23]. 

Nerve  stimulation  does  have  one  advantage.  The  probe  placement 
can  be  checked  by  observing  the  patient's  reaction  to  the  stimulating 
effect  of  the  current.  However,  other  nerves,  as  well  as  those  in  the 
lesion,  could  also  be  stimulated. 

Herrero  [14]  states  another  disadvantage  of  the  dc  method: 
"Animals  with  dc  lesions  usually  had  extensive  scar  tissue  both  around 
and  infiltrated  with  the  tissue". 

It  should  also  be  noted  that  the  dc  method  depends  on  both  the 
thermal  and  electrical  properties  of  the  tissue.  The  thermal  properties 
are  known  within  an  accuracy  of  approximately  5%  [7]  and  are  relatively 
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constant  throughout  the  tissue.  However,  the  electrical  properties  are 
not  well  known,  are  anisotropic  and  tend  to  vary  with  location. 

One  advantage  of  the  dc  method  is  that  the  electrodes  are  small; 
approximately  two  millimeters  in  diameter.  Therefore  the  probes  can 
be  painlessly  inserted  and  removed,  and  cause  little  damage  to  the 
surrounding  tissue. 

In  summary,  although  the  dc  method  can  produce  localized  lesions, 
the  irregularity  of  the  lesions,  control  difficulties,  and  possible 
side  effects  of  the  current,  greatly  diminish  its  usefulness. 

B.   RADIO-FREQUENCY  (rf)  CURRENT  METHOD 

The  radio-frequency  technique  produces  lesions  in  much  the  same  way 
as  the  direct  current  method.  A  long,  thin,  needle-like  probe  emits 
low  power  continuous  wave  radio-frequency  current.  An  indifferent 
electrode  is  placed  on  some  other  part  of  the  body  and  acts  as  a  ground. 
As  current  is  passed  between  the  electrodes,  the  resistivity  of  the 
tissue  causes  electrical  energy  to  be  converted  to  tissue  internal 
energy,  thus  increasing  the  local  tissue  temperature.  As  in  the  dc 
method,  lesions  are  a  function  of  the  tissue's  thermal  and  electrical 
properties. 

Gengler  [10]  presents  a  detailed  analysis  of  the  rf  technique. 
In  his  summary  of  the  method  he  states:  "The  radio-frequency  method 
allows  controlled  and  predictable  lesions  to  be  produced  within  safe 
limits  for  physiological  and  clinical  applications." 

Unlike  direct  current,  rf  current  has  no  preferred  path  in  tissue 
and  no  stimulating  effect.  However,  low-frequency  alternating  current 
does  have  a  stimulating  effect  in  tissue.  Therefore,  the  probe  location 
can  be  checked  by  passing  low-frequency  current  between  the  electrodes 
and  observing  the  patient's  response. 
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Brodkey  and  others  [3]  discuss  another  way  of  assuring  proper 
probe  location.  They  found  that  a  temperature  rangeof  40-49°C  reversibly 
blocks  nervous  function.  Thus,  after  positioning  the  probe,  and  raising 
the  local  tissue  temperature  to  the  reversible  range,  the  patient's 
response  can  indicate  if  the  probe  is  positioned  correctly. 

Like  the  dc  electrode,  the  radio-frequency  probe  is  small,  can  be 
painlessly  inserted  and  removed,  and  causes  little  damage  to  the 
surrounding  tissue. 

The  rf  method  does  have  several  disadvantages.  The  local  tissue 
temperature  may  exceed  100°C  causing  boiling, gas  formation,  carboniza- 
tion of  the  tissue,  and  adherence  of  coaguluum  to  the  probe  [2].  Also, 
defects  in  the  probe  insulation  have  been  known  to  cause  irregular 
lesions  along  the  probe  tract.  However,  these  difficulties  can  be 
avoided  by  carefully  controlling  the  various  physical  parameters  used 
in  the  rf  technique. 

In  summary,  the  radio-frequency  method  can  be  used  to  produce  well - 
circumscribed  lesions  and  the  probe  location  can  be  checked  by  two 
methods.  However,  as  in  the  dc  method,  lesions  are  still  a  function  of 
both  the  thermal  and  electrical  properties  of  the  tissue. 

C.   HIGH  RESISTANCE  HEATING  PROBE  METHOD 

In  this  method,  a  high  resistance  electrode  is  inserted  into  the 
tissue.  Current  is  passed  through  the  electrode  and  the  resistance  of 
the  electrode  causes  the  probe  temperature  to  increase.  Heat  is 
conducted  away  from  the  probe  into  the  surrounding  tissue.  A  region 
of  destroyed  tissue  is  formed  when  the  local  tissue  temperature  is 
raised  above  approximately  55°C  [2]. 
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Theoretically,  there  is  less  uncertainty  in  producing  lesions 
purely  by  heating  than  there  is  in  the  direct  current  or  radio-frequency 
techniques.  The  dc  method  produces  both  thermal  and  electrolytic  tissue 
destruction.  As  the  lesion  forms,  it  offers  an  increased  resistance  to 
the  passage  of  current.  Thus,  the  parameters  are  changing  as  the  lesion 
forms.  The  rf  method  minimizes  electrolysis  but  still  depends  on  the 
electrical  and  thermal  properties  of  the  tissue.  On  the  other  hand,  the 
high  resistance  probe  depends  only  on  the  tissue's  thermal  properties. 

Carpenter  and  Whittier  [4]  attempted  to  produce  lesions  in  animals 
by  using  a  heating  electrode.  Although  their  results  were  poor,  they 
concluded  that  it  was  due  to  the  crude  instruments  employed.  They 
stated:  "It  would  seem  profitable  to  explore  this  technique  with  a  more 
finely  perfected  thermocautery." 

Gildenberg  [11]  conducted  studies  with  a  high  resistance  heating 
electrode  and  a  rf  probe  and  found  that  the  heating  probe  produced  less 
variable  lesions  than  the  rf  probe.  In  addition,  the  rf  lesions  had 
significantly  more  lesions  along  the  puncture  canal  than  lesions  formed 
by  the  heating  probe. 

As  discussed  earlier,  a  temperature  range  of  40-49°C  can  reversibly 
block  nervous  function.  This  technique  can  be  used  to  check  the 
position  of  the  heating  electrode. 

Since  lesions  formed  using  the  resistance  heating  method  depend 
only  on  the  temperature  field  surrounding  the  probe,  the  surface 
temperature  of  the  probe  must  be  accurately  measured  and  controlled. 
This  has  been  orfe  of  the  difficulties  of  previous  resistance  probe 
designs.  Watkins  [23]  used  a  probe  with  a  thermocouple  placed  in  the 
tip.  He  found  that  the  actual  temperature  of  the  probe  could  exceed 
the  recorded  temperature  by  several  degrees  if  the  thermocouple  was  not 
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positioned  carefully.  Thermocouples  placed  in  the  surrounding  tissue 
to  record  the  temperature  can  distort  the  temperature  field  and  give 
erroneous  readings.  Gildenberg  [11]  inserted  a  thermocouple  one 
millimeter  from  the  heating  electrode  to  record  the  tissue  temperature. 
He  determined  that  controlling  the  current  to  maintain  a  constant 
temperature  at  the  thermocouple  was  no  more  accurate  than  maintaining 
a  constant  current. 

It  should  also  be  noted,  that  since  the  lesions  produced  by  the 
heating  probe  depend  only  on  thermal  destruction  of  the  tissue,  the 
lesions  are  smaller  than  those  produced  by  the  direct  current  or  radio- 
frequency  current  method. 

In  summary,  a  properly  designed  high  resistance  heating  probe  can 
produce  small,  safe  and  predictable  lesions.  The  uncertainty  is  less 
than  the  radio-frequency  or  direct  current  method  since  tissue  destruction 
depends  only  on  the  thermal  properties  of  the  tissue. 
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III.  HIGH  RESISTANCE  HEATING  PROBE 

A.   THEORY 

It  is  well  known  that  the  flow  of  electric  current  against  any 
resistance  is  accompanied  by  a  dissipation  of  electrical  energy  in  the 
form  of  heat.  The  amount  of  heat  generated  is  expressed  by  Joule's 
law: 

Q  =  IV  =  I2R  (1) 

Q  =  heating  rate 

V  =  voltage 

I  =  current 

R  =  resistance 
The  resistance  of  a  conductor  at  ordinary  temperature  is  a  function 
of  the  length  and  the  cross-section  of  the  conductor,  varying  directly 
with  the  length  and  inversely  with  the  cross-section  according  to  the 
following  relationship: 

R=P'f  (2) 

R  =  resistance  observed 
l  -   length  of  sample 
A  =  cross-section  area  of  sample 
The  proportionality  constant  p\   is  called  the  resistivity  of  the 
material.  Therefore,  all  other  factors  remaining  equal,  the  material 
having  the  highest  resistivity  will  generate  the  most  heat. 

All  metals  exhibit  some  change  in  resistance  with  temperature. 
For  pure  metals,  there  is  a  definite  and  constant  relationship  between 
the  temperature  and  the  resistance  of  the  metal.  This  characteristic 
is  used  in  electrical -resistance  thermometers  where  temperature  is 
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indicated  by  measuring  the  change  of  resistance.  Since  electrical 
resistance  measurements  can  be  made  with  high  accuracy  on  bridge 
circuits,  this  type  of  instrument  can  give  precise  temperature 
measurements. 

The  theory  of  the  resistance  thermometer  was  applied  to  the  design 
of  the  high  resistance  heating  probe.  From  the  above  discussion,  a 
given  metal  will  have  a  particular  resistance  at  a  particular  temperature 
The  inverse  of  this  was  used  in  designing  the  probe.  That  is,  imposing 
a  particular  resistance  on  a  metal  should  cause  the  metal  to  come  to  the 
corresponding  temperature  for  that  resistance.  The  resistance  can 
easily  be  controlled  by  regulating  the  voltage  and  current.  Thus,  the 
temperature  of  the  probe  can  be  measured  and  controlled  accurately 
without  the  use  of  any  external  measuring  devices. 

The  change  in  resistance  with  temperature  of  a  metal  is  given  by 
the  following  expression: 

R2  =  R1  [1  +  a'  (T2  -  T^]  (3) 

R-,  =  resistance  at  temperature  T, 
R«  =  resistance  at  temperature  T~. 

The  constant  a',  is  the  temperature  coefficient  of  resistivity  and 
is  expressed  as  ohms  per  ohm  per  degree.  It  should  be  noted  that,  in 
general,  a'  is  a  function  of  temperature. 

B.   DESIGN 

Equations  (1),  (2)  and  (3)  express  several  factors  that  were  con- 
sidered in  designing  the  high  resistance  probe.  From  equations  (1) 
and  (2)  it  can  be  seen  that  resistivity  should  be  high  to  generate  the 
most  heat.  Equation  (3)  shows  that  the  temperature  coefficient  of 
resistivity  should  also  be  high  so  that  the  change  in  resistance  with 
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temperature  will  be  large  and  therefore  easier  to  measure  and  control. 
Ideally,  a'  should  be  constant  over  the  range  of  temperatures  to  be 
measured. 

Several  metals  were  investigated.  Their  properties  are  listed  in 
Table  I.  In  addition,  the  properties  of  thermistors  were  also 
investigated.  Although  they  are  small,  rugged,  can  be  made  in  any  shape 
and  have  a  long  life,  their  change  in  resistance  with  temperature  is 
extremely  nonlinear.  Since  this  makes  temperature  control  difficult, 
thermistors  were  not  considered  further. 

TABLE  I 


Metal 

Resistivity 
(at  20°C) 

Coefficient 
of  Resistivity 

Constant 

Platinum 

10  microhm/cm 

.0038  ohm/ohm/°C 

Nickel 

7.8 

.006 

Ni chrome 

100 

.0004 

Copper 

1.72 

.0038 

Aluminum 

2.8 

.0039 

Balco 


Iron 


19.9 


10 


0045 


005 


25-100°C 


After  comparing  the  properties  of  the  metals  in  Table  I,  Balco  was 
chosen.  Balco  is  the  commercial  name  of  an  alloy,  consisting  of  70% 
nickel  and  30%  iron,  that  has  the  properties  that  are  desirable  in  a 
heating  electrode.  Its  temperature  coefficient  of  resistivity  is  large 
(a'  =  19.9  microhm/cm)  and  nearly  constant  between  25°C  and  100°C  [24]. 

Electrode  insulation  was  an  important  consideration  in  the  probe 
design,  since  a  voltage  leak  in  the  probe  can  cause  tissue  electrolysis 
Balco  is  available  with  a  standard  enamel  insulation.  The  enamel  is 
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made  from  an  oil  modified  expoxy  based  resin  and  has  a  continuous  duty 
operating  temperature  of  150°C  [24]. 

Various  shaped  heating  probes  were  considered.  However,  to 
facilitate  the  construction  of  the  probe,  a  cylindrical  shape  was  chosen. 
A  probe  length  to  diameter  ratio  (L/D)  of  ten  to  one  was  used  to  reduce 
end  effects  and  to  ensure  good  experimental  data. 

In  order  to  increase  the  extent  of  the  temperature  field  in  the 
test  medium,  thus  increasing  the  accuracy  of  the  experimental  data,  the 
probe  was  made  approximately  four  times  larger  than  an  actual  surgical 
probe. 

The  heating  electrode  (see  Figures  1  and  2)  consists  of  a  skeleton 
of  0.3894  centimeter  (cm)  diameter  glass  tubing.  Two  0.038  cm  copper 
wires  were  inserted  in  the  tubing.  One  wire  was  passed  out  the  tip 
and  high  resistance  0.0068  cm  Balco  wire  was  attached.  The  Balco  wire 
was  tightly  wound  back  along  the  glass  tubing  for  approximately  3.9  cm. 
It  was  then  led  back  into  the  tubing  through  a  small  hole  and  attached 
to  the  second  copper  wire.  The  holes  at  the  tip  and  side  of  the  glass 
tubing  were  closed  with  expoxy.  The  probe  was  then  given  a  thin 
coating  of  General  Electric  No.  7031  Adhesive  and  Insultating  Varnish 
thinned  with  Xylene  to  ensure  against  any  insulation  defects  in  the 
Balco  wire.  The  outside  diameter  of  the  completed  probe  was  0.4145  cm. 

C.   TEMPERATURE  CONTROL 

Precision  resistance  measurements  can  be  made  by  using  a  Wheatstone 
bridge.  The  bridge  (see  Figure  3)  consists  of  four  resistance  arms,  a 
voltage  source  and  a  detector.  When  A  is  to  B,  as  X  is  to  the  standard, 
there  will  be  no  potential  difference  bewteen  points  d  and  b.  Therefore, 
whenever  the  detector  reads  zero  voltage  from  point  d  to  b,  the  bridge 
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Figure  1.  Construction  of  Probe 
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Figure  3.  Wheatstone  Bridge. 
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is  in  balance  and  the  unknown  resistance  (X)  is  obtained  by  using  the 
relationship: 

X  =  (  |  )  x  (Standard)  (4) 

The  probe  was  calibrated  using  a  Rosemont  variable  temperature  oil 
bath  and  commutating  bridge.  The  calibration  curve  is  shown  in  Figure 
5.  Because  of  the  experimental  apparatus,  a  probe  temperature  between 
45°C  and  50°C  was  desired.  The  reason  for  this  choice  will  be  discussed 
later.  It  should  be  noted,  however,  that  in  actual  use,  a  probe 
temperature  between  85°C  to  95°C  would  be  used. 

As  discussed  earlier,  the  flow  of  current  against  a  resistance 
causes  heat  to  be  generated.  There  is  also  a  change  in  resistance  with 
temperature.  Since  the  temperature  of  the  probe  must  be  controlled 
precisely,  it  was  necessary  to  ensure  that  the  bridge  resistances  would 
not  change  as  heat  was  generated  in  them.  Therefore,  56.2  ohm  resistors 
rated  at  five  watts  were  used  for  resistances  A  and  B.  Referring  to 
equation  (4),  it  can  be  seen  that  since  the  ratio  of  A  to  B  is  one,  the 
standard  resistance  will  also  be  the  resistance  of  the  unknown  (probe). 

The  standard  resistance  was  wound  using  Evanohm  wire.  Evanohm  is 
a  commercial  wire  used  in  precision  resistance  standards  and  is  an  alloy 
of  75%  nickel,  20%  chromium,  2.5%  aluminum  and  2.5%  copper.  It  has  a 
low  temperature  coefficient  of  resistance  over  a  wide  range  of 
temperatures  [24].  The  completed  Evanohm  standard  had  a  resistance  of 
356.33  ohms  which  corresponded  to  a  probe  temperature  of  48.5°C  (see 
Figure  5). 

The  final  probe  circuit  (see  Figure  4)  consisted  of  the  two  56.2 
ohm  resistors,  the  Evanohm  standard  and  the  probe.  A  Lambda  Model  LP 
413  FM  regulated  power  supply  was  used  to  provide  power  to  the  probe. 
A  Weston  Model  1240  digital  voltmeter  was  used  as  the  detector. 
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Since  the  probe  had  a  high  temperature  coefficient  of  resistance, 
the  probe  resistance  was  constantly  changing  as  it  generated  heat. 
However,  by  manually  regulating  the  Lambda  power  supply  the  bridge  could 
be  kept  balanced  and  a  probe  temperature  of  48.5°C  maintained. 
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IV.  EXPERIMENTAL  APPARATUS 

A.   TEST  MEDIUM 

A  region  of  nervous  tissue  is  quite  complex.  However,  to  test  the 
heating  probe,  a  simplified  medium  that  had  the  heat  transfer  character- 
istics of  tissue  was  desired.  Since  the  lesions  produced  by  the  heating 
probe  are  a  result  of  thermal  conduction,  heat  transfer  by  convection 
had  to  be  negligible  in  the  test  medium.  In  addition,  the  medium  had 
to  be  clear  so  that  the  region  surrounding  the  probe  could  be  seen. 

Grayson  [13]  suggested  that  the  water  content  is  a  major  factor 
in  determining  the  thermal  conduction  of  animal  tissue.  Spells  [22] 
investigated  thermal  conduction  in  both  fluids  and  tissue  and  found 
there  was  a  linear  relationship  between  the  conductivity  and  the  water 
content  of  the  tissue.  The  linear  relationship  was  particularly  close 
when  the  water  content  was  above  50%. 

Cooper  and  Trezek  [7]  obtained  results  that  indicated  the  thermal 
properties  of  human  tissue  can  be  predicted  within  5%  by  measuring  the 
water  content  and  assuming  the  remaining  portion  of  the  tissue  is 
composed  of  equal  amounts  of  fat  and  protein.  They  also  found  that  the 
percentage  of  water  in  human  tissue  varied  between  71%  and  84%. 

The  above  results  indicate  that  an  ideal  test  medium  should  have  a 

high  percentage  of  water.  Gels  would  appear  to  be  perfect.  Gel 

consisting  of  gelatin  mixed  with  an  appropriate  percentage  of  water 

would  approximate  tissue  as  a  mixture  of  protein  and  water,  and  thermal 

properties  close  to  those  of  human  tissue  could  be  expected.  However, 

a  characteristic  of  gels  is  that  they  melt  at  approximately  22°C.  Since 

this  would  lead  to  convection  effects,  gel  was  not  used  as  the  test 

medium. 

29 


It  was  found,  however,  that  agar-water  solutions  consisting  of  as 
little  as  0.25%  agar  are  clear,  solidify  at  room  temperature,  and  remain 
solid  up  to  a  temperature  of  approximately  85°C.   In  addition,  the  thermal 
properties  of  a  1%,  or  less,  agar-water  solution  are  almost  identical 
with  those  of  pure  water,  and  within  10%  of  the  thermal  properties  of 
most  human  tissues  [7].  For  these  reasons,  a  mixture  of  97.7%  water  and 
0.3%  Difco  Bacto-Agar  was  chosen  as  the  experimental  test  medium. 

B.   LIQUID  CRYSTALS  AND  TEMPERATURE  MEASUREMENT 

Under  appropriate  conditions  of  temperature  and  pressure,  matter 
can  exist  in  a  solid  state,  isotropic  liquid  state  or  gas  state.  The 
majority  of  organic  and  inorganic  compounds  are  comprised  of  a  regular, 
three-dimensional  array  of  molecules  or  ions  in  the  solid  state,  i.e., 
a  crystal . 

When  the  crystal  lattice  is  heated,  thermal  agitation  makes  it 
impossible  for  the  molecules  to  maintain  a  regular  arrangement.  At  a 
given  temperature,  the  crystal  begins  to  melt  and  changes  to  a  liquid. 
This  change,  known  as  the  melting  point,  is  usually  very  sharp  and  a 
characteristic  of  the  compund.  At  still  higher  temperatures,  the  liquid 
boils  and  passes  into  the  gas  phase. 

Near  the  turn  of  the  century  it  was  discovered  that  some  compounds 
had  a  well  defined  melting  point  but  the  melt  was  opaque.  The  opacity 
disappeared  when  the  temperature  was  increased  to  a  definite  limit, 
given  a  true  isotropic  liquid.  Although  this  intermediate  phase  was 
fluid  to  some  extent,  it  also  exhibited  some  characteristics  of  a 
crystalline  solid.  This  phase  was  called  a  mesophase  and  compounds 
exhibiting  this  phenomenon  became  known  as  liquid  crystals  [12]. 
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There  are  three  types  of  mesophase  [12]. 
(i)  smectic  -  molecules  have  a  parallel  orientation  and  are  in  layers 
with  their  axes  at  right  angles  to  the  layers. 

(ii)  namatic  -  molecules  have  a  parallel  orientation  of  their  axes, 
(iii)  cholesteric  -  molecules  have  a  parallel  orientation  and  are  in 
layers  with  their  axes  in  the  plane  of  the  layer. 

The  cholesteric  mesophase  is  quite  different  from  the  smectic  and 
nematic  and  possesses  unique  optical  characteristics.  Cholesteric 
liquid  crystals  having  a  certain  structure  exhibit  dramatic  changes  in 
color  with  \/ery   small  changes  in  temperature.  This  color-temperature 
relationship  is  constant  and  reappears  everytime  the  liquid  crystal 
passes  through  a  particular  temperature  range. 

Dowden  [9]  discusses  the  structure  of  cholesteric  liquid  crystals. 
These  crystals  are  derived  from  cholesterol  or  other  sterol  systems 
and  have  a  basically  planar  structure.  The  axial  direction  in  each 
successive  layer  is  displaced  slightly  from  that  in  the  previous  layer 

o 

giving  a  spiral  trace.  The  thickness  of  each  layer  is  about  3A  and 

o 

the  pitch  of  the  helix  is  about  5500A  . 

The  property  which  is  the  basis  of  color  behavior  is  circular 
dichroism.  An  incident  beam  of  unpolarized  light  is  split  into  two 
components  having  electrical  vectors  rotated  in  opposite  directions. 
One  component  is  transmitted  and  the  other  is  scattered.  The  scattered 
light  normally  has  a  wavelength  peak  that  occurs  in  the  visible  region 
of  the  spectrum. 

Temperature  changes  disturb  the  intermolecular  forces  and  cause  a 
shift  in  the  molecular  structure.  This  changes  the  wavelength  of  the 
light  that  is  scattered.  Thus,  cholesteric  liquid  crystals  exhibit 
dramatic  changes  in  color  with  very  small  changes  in  temperature. 
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Since  the  colors  scattered  by  liquid  crystals  are  unique  for  specific 
temperatures,  the  quantitative  measurement  of  temperature  is  possible 
to  an  accuracy  of  less  than  0.1 °C  [8]. 

Liquid  crystals  can  be  applied  directly  to  a  surface.  However,  a 
dark  uniform  layer  is  normally  placed  on  the  surface  before  the  liquid 
crystals.  This  provides  a  standard  surface  texture,  intensifies  the 
resulting  colors  and  isolates  the  crystals  from  surface  contamination. 

Liquid  crystals  have  some  application  problems.  Since  they  are 
viscous,  applied  film  thicknesses  can  easily  vary  and  result  in  non- 
reproducible  thermal  patterns.  Direct  exposure  of  the  crystals  to 
solvent  vapors,  oils,  dirt  and  ultra-violet  light  can  cause  variations 
in  their  sensitivity.  However,  it  is  possible  to  encapsulate  the  liquid 
crystal  in  a  tiny  capsule  (10  to  30  microns),  thus  shielding  it  from 
most  of  these  detrimental  effects  [17]. 

Encapsulated  crystals  are  available  commercially  in  a  variety  of 
forms  that  cover  a  wide  range  of  temperatures.  Rochrome  liquid  crystal 
tapes  manufactured  by  Hoffman-LaRoche  Inc.  were  used.  These  crystals 
were  coated  on  a  mylar  subtratum  with  an  absorptive  black  background. 
The  overall  thickness  of  the  crystals  and  mylar  was  less  than  0.01 
inches. 

The  Rochrome  crystals  are  calibrated  with  an  accuracy  of  0.1 °C  and 
have  a  response  time  of  less  than  one  second  [18].  Parker  [19]  found 
that  the  response  time  of  the  encapsulated  crystals  depends  on  the 
crystal  thickness.  For  a  thickness  of  0.0035  inches,  response  time  was 
of  the  order  of  200  to  300  milliseconds. 

By  placing  the  liquid  crystal  tape  in  the  test  medium,  the  tempera- 
ture field  around  the  probe  could  be  seen  and  temperatures  throughout 
the  field  accurately  measured. 
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Rochrome  liquid  crystal  tape  lot  no.  AS-70-36  was  used  in  the  actual 
experiments.  The  color  to  temperature  relationship  of  this  tape  was: 

29.9°C  -  red 

31.5°C  -  green 

34.2°C  -  blue 
During  preliminary  experiments  with  the  liquid  crystal  tapes  it  was 
found  that  the  water  in  the  test  medium  caused  the  encapsulated  crystals 
to  separate  from  the  mylar  backing.  To  prevent  this  from  happening,  the 
tapes  were  given  a  thin  coating  of  Varathane  and  the  edges  were  sealed 
with  a  rubber  adhesive. 

C.   TEST  CELLS 

Two  test  cells  were  used  to  study  the  three-dimensional  temperature 
field  surrounding  the  probe  (see  Figure  6).  A  circular  container,  five 
inches  in  diameter  and  four  inches  deep,  was  constructed  of  plexiglass 
and  was  used  to  visually  observe  the  radial  temperature  field.  A 
removable  ring  was  fitted  into  the  cell  to  support  the  liquid  crystal 
tape  (see  Figure  7).  The  tape  was  attached  to  a  thin  sheet  of  mylar 
that  was  streched  across  the  ring.  The  probe  was  inserted  through  the 
bottom  of  the  cell  and  was  held  in  place  by  a  bushing  that  allowed  the 
height  of  the  probe  to  be  adjusted. 

The  second  cell,  a  4  5/8  x  5  1/2  x  6  inch  rectangular  container, 
also  constructed  of  plexiglass,  was  used  to  observe  the  axial  temperature 
field  near  the  probe.  The  probe  was  inserted  through  the  side  of  the 
container.  A  removable  bracket  was  manufactured  to  hold  the  liquid 
crystal  tape  lengthwise  along  the  probe  (see  Figure  7).  Because  of 
symmetry,  the  temperature  field  on  only  one  side  of  the  probe  was 
measured. 
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To  ensure  that  the  container  walls  would  not  disturb  the  tempera- 
ture field  around  the  probe,  both  containers  were  constructed  so  that 
their  dimensions  were  large  compared  to  the  probe  dimensions. 
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V.  EXPERIMENTAL  PROCEDURE 

The  test  cells  were  filled  with  the  0.30%  agar-97.79%  water  test 
medium  for  the  experimental  runs.  The  probe  was  positioned  in  the 
circular  test  cell  so  that  there  was  an  equal  length  of  probe  above  and 
below  the  liquid  crystal  tape.  Thus,  the  temperature  field  around  the 
midpoint  of  the  probe  was  measured.  The  probe  was  covered  by  the  agar- 
water  medium  so  that  the  top  of  the  probe  was  at  least  five  probe 
diameters  below  the  surface. 

All  test  runs  were  conducted  at  room  temperature  using  Rochrome 
liquid  crystal  tape  lot  no.  AS-70-36  to  measure  the  temperature  field 
around  the  probe.  As  mentioned  earlier,  this  tape  had  the  following 
temperature  to  color  relationship: 

29.9°C-  -  red 

31.5°C  -  green 

34.2°C  -  blue 
A  narrow  strip  of  Series  9  Rochrome  tape  was  placed  in  the  test 
medium  to  measure  the  medium's  ambient  temperature.  The  temperature  to 
color  relationship  of  the  Series  9  tape  was: 

23.0°C  -  red 

24.0°C  -  green 

27.0°C  -  blue 
The  ambient  room  temperature,  normally  about  23.5°C,  was  used  as  a 
rough  check  of  the  temperature  obtained  from  the  Series  9  tape. 

Before  the  containers  were  filled  with  the  agar-water  test  medium, 
four  reference  marks,  4.57  centimeters  apart,  were  placed  on  the 
corners  of  the  tape,  Figure  8  is  a  top  view  of  the  circular  test 
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cell  showing  the  probe  in  place,  the  narrow  strip  of  Series  9  tape, 
and  the  four  reference  marks. 

The  resistance  of  the  probe  was  controlled  manually  by  adjusting 
the  Lambda  power  supply  to  maintain  zero  voltage  on  the  detector. 
Figure  9  is  an  overall  view  of  the  apparatus. 

Photographs  of  the  temperature  field  surrounding  the  probe  were 
taken  using  a  Graflex  4x5  Press  camera  fitted  with  a  Polaroid  Land 
film  holder.  Polaroid  4x5  Land  film  Type  57,  3000  speed  was  used. 

Photographs  were  taken  as  various  increments  of  time.  The  outer 
radius  of  the  temperature  field,  corresponding  to  29.9°C,  and  the 
reference  marks  were  measured  in  each  picture.  Using  the  actual 
distance  between  the  reference  marks  (4.57  cm)  and  the  probe  radius 
(0.2076  cm),  the  non-dimensional  radius  was  calculated. 
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VI.  THEORETICAL  ANALYSIS 

A.   ONE-DIMENSIONAL  MODEL 

1 .  Metabolic  and  Blood  Flow  Effects  Neglected 

In  order  to  simulate  the  heat  transfer  process  in  the  region 
surrounding  the  probe,  a  one-dimensional  heat  transfer  model  must  be 
developed.  In  this  model  it  will  be  assumed  that  the  tissue  is  an 
infinite  homogenous  medium  with  known  conductivity,  density  and  specific 
heat  and  that  it  is  initially  at  a  uniform  temperature  of  T  .  The 
probe  is  assumed  to  be  a  cylinder  of  infinite  length  and  radius  r  with 
a  constant  temperature  of  T  (see  Figure  10).  Thus,  the  problem  is  to 
solve  for  the  temperature  distribution  in  a  region  bounded  by  a  circular 
cylinder. 

The  general  conduction  equation  is:. 

£  Ar     fr)  =  pc|I  (5) 

k  =  thermal  conductivity 

p  =  density 

c  =  specific  heat 

a  =  —  =  thermal  diffusivity 
pc 

It  is  now  convenient  to  introduce  the  following  nondimensional 

variables: 

R  =  —  (distance) 
ro 

at   / .  •   \ 
t  =  ~2~   (time) 

ro 
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Figure  10.  Mathematical  Model  of  Probe 
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T  -  T 

e  =  j j-    (temperature) 

P  "  o 


Substituting  these  variables  into  equation  (5)  yields, 
1  3(R  air)    ae 


\3K/_  i  ~\ 


R    dR  3t 

The  boundary  conditions  are: 

(a)  at  the  probe  surface  (r  =  a,  or  R  =  1 ) ,  T  =  T  and  e  =  1. 

(b)  at  large  distances  from  the  probe  the  temperature  of  the  tissue 

remains  constant;  that  is ,  as  r  ->  °°,  T  -*■  T  or  nondimensionally, 

as  R  -*■  <»,  e  ->  o. 

The  initial  condition  is  at  time  t=o,T=T  or  e  =  0  for  all 

o 

locations. 

Carslaw  and  Jaeger  [5]  give  the  closed  form  solution  of  equation 
(6): 

e  .  i  .  1  f"   e  -T"2C0  (u.Ru)du 

tt  nJ = (7) 

°    U  [J2  (u)  +  Y2  (u)] 

CQ  (u,  Ru)  =  JQ  (u)  YQ  (Ru)  -  Yq  (u)  Jq  (Ru) 

where,  J  and  Y  are  Bessel  functions.  Equation  (7)  was  evaluated 

numerically  by  Jaeger  for  various  values  of  t  and  R  [15]. 

2.  Metabolism  and  Blood  Flow  Effects  Included 

Pennes  [20]  was  the  first  to  account  for  metabolic  and  blood 

flow  effects  in  the  heat  equation.  Pennes  suggests  that  the  equation 

take  the  form: 

3/r  3T\ 
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where 


m 


.  =  blood  flow  rate  per  unit  volume 


S   =  metabolic  heat  rate  per  unit  volume 
T.  =  blood  temperature 

In  this  model,  the  effect  of  blood  flow  is  treated  as  a  perfect 
heat  exchanger.  That  is,  it  is  assumed  that  the  blood  enters  a  unit 
volume  of  tissue  at  the  arterial  temperature,  T.  ,  and  acquires  the 
tissue  temperature  as  it  exits  the  volume. 

If  the  metabolic  heat  generation  rate  is  assumed  to  be  constant,  it 
can  be  combined  with  the  blood  flow  term  as: 

mb  cb  <To  ~T) 

where 

S 

T  =  TL  +    m 


o    b    mb  cb 

This  corresponds  to  redefining  the  reference  temperature.  Cooper 
and  Trezek  [6]  show  that  for  brain  tissue,  T  is  on  the  order  of  0.5°C 
higher  than  the  arterial  blood  temperature,  T.  . 

Equation  (8)  may  be  normalized  by  defining  a  nondimensional  "blood 

flow  parameter": 

m  c  r 
@  =   b  r  °   (blood  flow  term) 

Substituting  this  variable  along  with  R,  t  and  e  into  equation  (8) 
yields: 


3/R  39\ 
1  \       3x/ 


6  6   =  !!-  (9) 


R    9R     p  u     3t 

The  boundary  conditions  and  initial  condition  are  the  same  as  were 
used  for  equation  (6). 

Equation  (9)  was  solved  with  the  use  of  Laplace  Transformations  and 
the  following  closed  form  solution  was  obtained: 
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2 
KIM)      o  _Rt  /-co  e  "TU  C  (u,  Ru)  du 

e=_^ ^_2eBr  r    o (10) 

Ko(/F)   *    o^   [u+e/u]  [J*  (u)  +  Y2  (u)] 

CQ  (u,  Ru)  =  JQ  (u)  YQ  (uR)  -  JQ  (uR)  YQ  (u) 

Equation  (10)  was  solved  numerically  for  various  values  of  R  and  t 
and  for  values  of  6  =  0.0,  0.1,  0.5  and  1.0  and  the  results  are  tabulated 
in  Tables  II  -  XI.  The  reader  is  referred  to  Appendix  A  for  details. 

B.   TWO-DIMENSIONAL  MODEL 

The  one-dimensional  model  assumed  that  the  probe  was  infinite  in 
length.  However,  the  probe  consists  of  a  constant  temperature  cylinder 
of  finite  length  and  a  cylindrical  insulated  stem.  The  one-dimensional 
solution  does  not  take  end  effects  into  account.  Thus,  it  was  desired 
to  have  a  two-dimensional  solution  that  would  include  these  effects. 

A  digital  computer  program,  by  the  code  name  of  Trump  [16],  was 
used  to  determine  the  two-dimensional  transient  temperature  field 
around  the  probe. 

Since  TRUMP  obtains  a  solution  in  dimensional  form,  a  nodal  netword 
(see  Figure  11)  had  to  be  constructed  for  the  geometric  configuration 
of  the  probe.  Due  to  symmetry,  it  was  necessary  to  generate  the  field 
for  only  one  half  of  the  cross-sectional  area.  The  dimensions  of  the 
network  were  chosen  large  enough  so  that  the  boundary  conditions  at 
R  =  30  and  X  =  24.0  r  did  not  affect  the  temperature  field  during  the 
transient  period  which  was  investigated.  The  network  therefore 
simulated  an  "infinite"  field. 

The  TRUMP  solution  gave  values  of  e  at  each  node  point  for  various 
values  of  nondimensional  time.  The  results  are  given  in  Tables  XII  - 
XVI. 
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Figure  11.     TRUMP  Nodal    Network. 


46    . 


VII.  RESULTS 

The  preceding  section  gave  details  of  the  three  models  that  were 
used;  a  one  and  two-dimensional  model  without  blood  flow  and  metabolic 
effects  and  a  one-dimensional  model  including  these  effects.  Numerical 
results  for  the  one-dimensional  models  are  given  in  Tables  II-XI. 
Tables  XII  -  XVI  give  the  numerical  results  for  the  two-dimensional 
model.  Values  for  e  are  given  at  various  nondimensional  locations  (R) 
and  times  (t).  The  experimental  results,  differences  between  the  one 
and  two-dimensional  models,  and  the  relative  effects  of  blood  flow  and 
time  will  now  be  discussed. 

A.   EXPERIMENTAL  RESULTS 

Figures  12,  13,  14,  15  and  16  are  representative  pictures  of  the 
radial  temperature  field  surrounding  the  probe  taken  at  t  =  0.0,  1.0, 
3.0,  5.0  and  10.0.  Figure  17,  18,  19,  20  and  21  depict  the  axial 
temperature  field  at  t  =  0.0,  1.0,  3.0,  5.0  and  10.0.  The  latter 
pictures  also  show  end  effects. 

It  was  convenient  to  use  the  following  nondimensional  variables  in 
calculating  the  experimental  results: 

T  -  T 
e  =  j — j-y—  (critical  temperature) 
P   o 

where,  T  is  the  temperature  corresponding  to  a  given  liquid  crystal 
color.  In  the  experimental  runs,  this  color  was  red  and  corresponded  to 
a  temperature  of  29.9°C.  In  a  surgical  application,  T  is  the 

temperature  at  which  tissue  destruction  begins. 

r 

R  =  —  (radius  of  lesion) 
c   rQ 

The  experimental  data  and  results  are  given  in  Table  XVII. 
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Figure  13.  Radial  Temperature  Field  for  t  =  1.0 
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Figure  14.  Radial  Temperature  Field  for  t  =  3.0, 
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Figure  15.  Radial  Temperature  Field  for  t  -  5.0, 
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Figure  16.  Radial  Temperature  Field  for  t  =  10.0, 
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Figure  17.  Axial  Temperature  Field  for  x  =  0.0. 
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Figure  18.  Axial  Temperature  Field  for  t  =  1.0, 
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Figure  19.  Axial  Temperature  Field  for  t  =  3.0, 
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Figure  20.  Axial  Temperature  Field  for  x  =  5.0, 
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Figure  21.  Axial  Temperature  Field  for  t  =  10.0, 
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The  experimental  results  were  compared  with  the  two-dimensional 
model.  Figure  22  is  a  typical  graph  of  the  nondimensional  time  for  a 
critical  temperature,  e  =  0.200.  The  solid  line  represents  the  two- 
dimensional  solution  and  the  dashed  lines  represent  the  error  bound. 
The  reader  is  referred  to  Appendix  B  for  details  of  the  error  analysis. 

The  experimental  results  fall  within  the  error  bounds  except  at 
short  times.  This  is  due  to  the  fact  that  the  analytical  solutions 
were  obtained  assuming  a  step  input.  That  is,  it  was  assumed  that  the 
probe  came  to  its  final  temperature  immediately.  Because  of  equipment 
limitations  this  was  impossible  to  duplicate  in  the  experimental  runs 
and  it  took  the  probe  approximately  15  seconds  to  reach  its  final 
temperature. 

Figure  23  is  a  graph  of  the  nondimensional  lesion  radius  versus 
nondimensional  time  for  a  critical  temperature  of  e  =  0.235.  The 
solid  line  represents  the  two-dimensional  solution. 

The  axial  temperature  field  is  shown  in  Figure  24.  Again,  the 
solid  lines  represent  the  two-dimensional  solution.  Only  one  experi- 
mental run  was  conducted  using  this  configuration.  Although  the 
experimental  isotherms  exhibit  the  same  shape  as  the  analytical  iso- 
therms, they  are  lagging  in  time.  A  great  deal  of  difficulty  was 
encountered  in  placing  the  liquid  crystal  sheet  in  a  uniform  plane  next 
to  the  probe  and,  in  fact,  the  edges  of  the  sheet  tended  to  sag  after 
a  short  period  in  the  agar-water  solution.  More  work  is  needed  on 
this  particular  aspect  of  the  experiment  since  experimental  errors  are 
obviously  large. 
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Figure  24 


Axial  Position  (x)  vs.  R  for  Various  Values  of 
t  and  er  =  0.285. 
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B.   ONE-DIMENSIONAL  MODEL 

1 .  Effect  of  Blood  Flow  on  the  Temperature  Field 

Figure  25  is  a  graph  of  the  nondimensional  temperature  versus 
the  radial  location  in  the  tissue  for  a  nondimensional  time  equal  to 
five  (corresponds  to  a  real  time  of  approximately  150  seconds  for  the 
experimental  probe).  Four  curves  are  presented:  the  dashed  line 
represents  no  blood  flow  effect  and  the  solid  lines  represent  b's  of 
0.1,  0.5  and  1.0.  Note  from  the  graph  that  blood  flow  tends  to  decrease 
the  temperature  field,  resulting  in  a  slightly  smaller  lesion.  To  gain 
a  feeling  for  the  magnitude  of  this  effect,  a  typical  value  of  8  will 
be  calculated.  A  typical  value  for  m.    c.  of  blood  flow  in  the  brain  is 
0.01  cal/cm  /sec/°C  [21]  and  a  typical  value  of  the  thermal  conductivity 
of  brain  tissue  is  0.0014  cal/cm/°C/sec  [7].  Therefore,  for  brain 
tissue,  8  is  approximately  equal  to  7  r^,  where  r  is  the  probe  radius 
expressed  in  centimeters.  The  experimental  heating  probe  had  a  radius  of 
approximately  0.2  cm;  thus,  8  would  be  of  the  order  of  0.28.  However, 
an  actual  probe  would  be  approximately  one  fourth  the  size  of  the  experi- 
mental probe  and  8  would  then  be  much  less  than  0.1.  Figure  26  shows 
that  values  for  blood  flow  in  this  range  decreased  the  temperature  only 
slightly. 

2.  Effect  of  Increased  Time  on  Temperature  Field 

Figure  26  is  a  graph  of  the  nondimensional  temperature  versus 
nondimensional  time.  One  can  see  that  at  short  times,  the  rate  of 
lesion  growth  is  large,  but  as  time  increases,  the  lesion  growth 
decreases  and  approaches  a  steady  size. 

Figure  27  illustrates  the  effect  of  blood  flow  at  various  times. 
In  this  graph  the  temperature  at  a  radius  of  2.5  is  plotted  against  time 
for  8  =  0.0,  0.1,  0.5  and  1.0.  At  short  times  the  blood  flow  effect 
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seems  to  be  negligible  but  as  time  is  increased,  the  effect  of  blood 
flow  becomes  more  noticeable.  Since  the  probe  radius  enters  the  blood 
flow  parameter  as  a  squared  term,  the  blood  flow  effect  would  become 
an  important  parameter  for  a  larger  probe  at  longer  times. 

C.   TWO-DIMENSIONAL  MODEL 

Results  for  the  two-dimensional  model  are  tabulated  in  Tables  XII  - 
XVI.  The  results  are  given  for  various  positions  along  the  probe 
length,  where  the  X=0  position  corresponds  to  the  midpoint  of  the  probe. 

The  two-dimensional  solution  was  used  to  compare  the  experimental 
results  with  the  theoretical  results.  Figure  24  shows  the  isotherms 
along  the  probe  length  for  various  nondimensional  times.  By  extending 
this  graph  it  can  be  seen  that  the  temperature  field  around  a  cylindrical 
probe  has  an  elliptical  shape.  The  dimensions  of  the  ellipsoid  can  be 
changed  by  varying  the  length  to  diameter  ratio  of  the  probe.  If  this 
ratio  is  made  small  enough  ,  the  shape  of  the  temperature  field  will 
approach  that  of  a  sphere. 

Figure  28  is  a  graph  of  the  axial  position  (X)  versus  the  radial 
location  in  the  tissue  for  a  nondimensional  time  equal  to  five  and  e  = 
0.300.  The  solid  lines  represent  the  two-dimensional  solution  for  an 
L/D  of  ten,  and  the  dashed  line  represents  the  one-dimensional  solution. 
Around  the  midpoint  of  the  probe,  the  two  solutions  are  in  close 
agreement.  However,  as  the  distance  from  the  mid-point  increases,  the 
one-dimensional  solution  is  no  longer  accurate. 
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Figure  28.     Axial    Position   (X)  vs  R  for  T  =  5.0 
and  e  =  0.300. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

Previous  experimental  methods  used  to  predict  lesion  size  and  shape 
have  generally  reduced  to  trial  and  error.  The  usual  procedure  has  been 
to  produce  lesions  in  animals,  vary  the  time  and  current,  and  analyze 
the  results.- 

The  theoretical  predictions  of  the  temperature  field  obtained  from 
the  two-dimensional  heat  transfer  model  of  the  probe  fall  well  within 
the  error  bounds  of  the  experimental  results.  Thus,  it  appears  possible 
to  analytically  predict  the  size  of  lesions  produced  with  the  constant 
temperature  resistance  probe. 

Experimental  isotherms  are  reproducible  and  exhibit  sharp  delimita- 
tion, particularly  at  short  times.  .The  probe  itself  is  safe,  easy  to 
construct  and  simple.  The  temperature  control  circuit  is  small, 
accurate  and  easy  to  operate. 

All  results  are  presented  in  nondimensional  form.  However,  it 
should  be  noted  that  both  the  experimental  and  theoretical  models  are 
simplified  and  the  results  have  been  verified  only  over  a  narrow 
temperature  range  and  only  for  a  cylindrical  probe  with  a  length  to 
diameter  ratio  of  ten. 

The  blood  flow  effect  is  found  to  be  an  important  parameter  in 
determining  the  temperature  field  in  tissue  if  a  large  probe  is  used. 
For  very  small  probes,  the  effect  is  marginal.  Including  the  effect  of 
blood  tends  to  slightly  decrease  the  temperature  throughout  the  field 
and  a  steady  state  or  equilibrium  is  reached  at  an  earlier  time  than  if 
blood  flow  is  neglected.  As  time  is  increased,  the  effect  of  blood 


68 


flow  becomes  more  predominant.  Thus,  if  heat  is  applied  for  any  sub- 
stantial period  of  time,  the  effect  of  blood  flow  is  felt  more  strongly. 

B.   RECOMMENDATIONS 

It  is  recommended  that  additional  experimental  work  be  carried  out 
with  the  heating  probe  over  higher  temperature  ranges  and  with  various 
probe  dimensions  and  L/D  ratios. 

As  mentioned  previously,  problems  were  encountered  in  attempting 
to  place  the  liquid  crystal  tape  uniformly  along  the  axis  of  the  probe. 
Accurate  axial  isotherms  were  therefore  not  obtained  in  the  present  set 
of  experiments.  Further  investigation  is  needed  on  this  aspect  of  the 
experiment. 

It  is  also  recommended  that  a  technique  for  automatically  controll- 
ing the  probe  surface  temperature  be  developed.  Two  highly  desirable 
features  of  such  a  control  system  would  be:  1)  an  ability  to  select  a 
range  of  probe  surface  temperatures,  and  2)  an  ability  to  pre-select 
the  duration  of  the  heating  period.  These  two  features,  together  with 
the  analytical  model,  will  allow  a  more  accurate  prediction  of  lesion 
size. 
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APPENDIX  A:  ONE-DIMENSIONAL  MODEL  SOLUTION 

Statement  of  Problem: 

A  cylindrical  probe  of  radius  r  is  embedded  in  tissue  which  is 
initially  at  a  uniform  temperature  T  .  Within  the  tissue  there  is  a 
volumetric  heat  generation  due  to  metabolism  of  strength  S  .  In 
addition,  there  is  blood  flow  in  the  medium  with  a  flowrate  m,  , 
specific  heat  c.  and  temperature  L.  Assume  that  the  length  to  dia- 
meter ratio  of  the  probe  is  large  enough  so  that  end  effects  may  be 
neglected  and  assume  the  tissue  is  homogeneous  and  infinite  in  extent. 
The  probe  is  at  a  constant  temperature  T  .  Also,  assume  that  at  an 
infinite  distance  from  the  probe,  no  temperature  change  occurs.  Solve 
for  the  tissue  temperature  assuming  constant  properties  and 
one  dimensional  radial  heat  flow. 

Notation: 


Tissue  thermal  diffusivity  cm  /sec 


Symbol  Description  Typical  Units 
a 

c  Tissue  specific  heat  cal/gm/°C 

c.  Blood  specific  heat  cal/gm/°C 

k  Thermal  conductivity  of  tissue  cal/cm/sec/°C 

b 


3 

mK       Volumetric  blood  flowrate  gm/cm  /sec 


r  Radial  location  in  tissue  cm 

r  Probe  radius  cm 
o 

3 

S  Metabolic  heat  generation  rate  cal/cm  /sec 
m 

t  Time  sec 

T  Tissue  temperature  °C 
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Symbol      Description  Typical  Units 

T.       Arterial  temperature  °C 

T        Initial  tissue  temperature  °C 

T        Probe  temperature  °C 

Nondimensional  Groups: 


T-Tc 


e  =  j — j-j-  (temperature) 


R  =  —  (distance) 
o 


t  =  —  (time) 

•I 

m  c  r 
3  =  b  b  °  (blood  flow) 
k 

S 

T„  =  T.  +  — -  (initial  temperature) 

o    d   m,  c. 


The  general  heat  conduction  equation  in  cylindrical  coordinates  is 
3/r  3T\ 

Substituting  the  above  nondimensional  quantities  into  equation 

(A.l)  yields: 

3/R  39 


3/K  39_\ 
1   \   3R/ 


39 


M=r  (A. 2) 


R   3R     M     3t 
Equation  (A. 2)  must  be  solved  in  order  to  describe  the  temperature 
of  the  tissue  as  a  function  of  time.  In  order  to  solve  the  problem, 
boundary  conditions  and  initial  conditions  must  be  stated.  They  are 
as  follows: 
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(a)  at  r  =  r  ,  (at  the  probe  surface),  T  =  T  or,  in  nondimensional 
form  at  R  =  1 ,  8  =  1 , 

(b)  as  r  ->  °°,  the  medium  remains  unaffected  as  r  -*•  »,  T  ■*■   T  or, 
in  nondimensional  form  as  R  ■*  °°,  e  -*■  o, 

(c)  at  t  =  o,  T  =  T  or,  in  nondimensional  form  at  t  =  o,  e  =  o. 
By  using  Laplace  Transformations,  a  closed  form  solution  to 

equation  (A. 2)  can  be  obtained: 

2 
K     /BR       „  e  "TU     C     (u,   Ru) 

e  -  — f rr-% —  du  (A-3) 

Ko/3- 


2     r     e  L0  (u>   Ru] 


C     =  J     (u)   Y     (uR)   -  J     (uR)   Y     (u) 

0  00  00 

If  the  blood  flow  term  is  neglected,  equation  (A. 3)  reduces  to: 

2 

r     e  "TU   c    (Uj  Ru)  ^ 

9  =  1-|   /    n °~ ~? (A.4) 

77  oJ        u  [Jj  (u)  +  Y^  (u)] 

Equation  (A.4)  is  the  closed  form  solution  given  by  Carslaw  and 
Jaeger  [5].  This  integral  has  been  solved  by  Jaeger  [15]. 


Solution: 

A  finite  increment  technique  was  used  to  solve  for  the  temperature 
field.     Equation   (A. 3)  was  rewritten  in  the  following  form: 


e  =  _o__) 1     _  £        /       q   (u)  du  (A. 5) 

K 


^J     -|       /      q    (u)   du 


e"TU     C     (u,  Ru) 
q  (u) 


(u  ♦  i)[jg   (u)   t  Y2   (u)] 

Assume  an  approximation  to  the  integral  can  be  derived  by  passing 
an  interpolating  function  through  several  adjacent  points  and  then 
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integrating  the  interpolating  function.  Simpson's  Rule  uses  the 

parabola  as  an  interpolatory  function  between  u._-,,  u.  and  u.  ,  where 

the  u-  are  evenly  spaced  values  of  the  interval    (u    ,  u   )  with  a 
l  on 

subinterval   length  h  =  u.+-,   -  u..     Integrating  the  parabolic 
interpolating  function  gives: 


(A. 6) 


ui+l 
u../  *»du4[,i-l   +4yi+Yi+l   3 

where  Y-   =  q(u. ) 

Expanding  equation  (A. 6)  over  the  entire  interval  yields, 

f      q(u)  du  -  5-  [Y0  +  4Y,  +  2Y0  +  . . .  4Y  ,  +  Y  ]         (A. 7) 
0J  3  L  o    1     2       n-1    nJ        K       ' 

Computer  programs  1,  2  and  3  were  used  to  solve  this  problem  for 
3  =  0.0.  Programs  4,  5,  and  6  were  used  to  solve  for  6  =  0.1.  To 
solve  for  other  values  of  B  only  one  card  must  be  changed. 

Because  q(u)  is  singular  at  u  =  0,  the  lower  limit  of  the  integral 
was  taken  as  0.001.  The  error  caused  by  this  was  assumed  to  be  negli- 
gible. By  varying  the  upper  limit  of  the  integral  between  1.0  and  100.0, 
it  was  found  that  the  solution  coverged  at  approximately  8.0.  Therefore, 
an  upper  limit  of  10.0  was  used. 

The  Y  Bessel  functions  involving  the  lower  limit  of  the  integral, 
0.001,  were  solved  using  a  separate  program  and  calling  the  BESY 
subroutine.  This  was  done  to  avoid  underflow  and  reduce  computer  time. 

The  included  programs  solve  the  temperature  field  for  times 
between  0.1  and  1.0  only.  To  change  the  time  interval  only  one  card 
must  be  changed. 

Jaeger's  results  [15]  were  used  as  a  check  on  the  accuracy  of  the 
results  obtained  from  Programs  1,  2  and  3  when  $  =  0.0.  The  greatest 
difference  between  the  solutions  was  0.001. 
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APPENDIX  B:  ERROR  ANALYSIS  OF  EXPERIMENTAL  RESULTS 


Three  nondimensional  quantities  were  calculated  in  analyzing  the 


experimental  results. 


Tc  -  T 

e     =  j y~    (critical   temperature) 

p  "  o 

r 
R„  =  —  (critical  radius) 
c   ro 

t  =  2Li  (time) 
ro 

Ambrosius  [1]  gives  the  following  equations  for  the  resultant 

errors  (R)  in  computations  involving  measurements: 

(a)  addition 

(A  ±  a)  +  (B  ±  b) 

R  =  ±  (a2  +  b2)2  (B.l) 

(b)  subtraction 

(A  ±  a)  =  (B  ±  b) 

R  =  ±  (a2  +  b2)2  (B.2) 

(c)  multiplication 
(A  ±  a)  (B  ±  b) 

R  =  ±  [(Ab)2  +  (Ba)2]2~  (B.3) 

(d)  division 

(B  ±  b)  /  (A  ±  a) 

.  ,  R  „  ±  r(Ba/A)2  +  (b)2  f  (B>4) 

where  A  and  B  are  two  measurements  and  an  error  of  ±  a  is  made  in  the 
measurement  of  A  and  an  error  of  ±  b  is  made  in  the  measurement  of  B 
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The  following  errors  were  estimated  for  the  variables  in  9  : 
3  c 

T  =  ±  0.1°C 
c 

TQ  =  ±  0.3°C 

T  =  ±  0.2°C 
P 

An  examination  of  e  along  with  equation  (B.4)  indicates  that  the 

largest  error  in  e  occurs  when  T  is  the  greatest,  or  when  the  quantity 

T  -  T  is  the  smallest.  The  highest  value  of  T  was  in  test  run  10. 
co  o 

The  resultant  errors  will  now  be  estimated  using  data  from  run  10. 

TQ  =  25.2°C  ±  0.3°C 

T  =  29.9°C  ±  0.1°C 
c 

T  =  48.5°C  ±  0.2°C 

T  -  T  =  4.7°C 
c    o 

Tn  -  "L  =   18.6°C 
P    o 

Substituting  these  values  into" equation  (B.2),  the  following  errors 

are  obtained: 

T  -  T  =  ±  0.32°C 
c   o 

T  -  T  +  ±  0.36°C 
P    o 

From  equation  (B.4)  the  resultant  error  in  9  is: 
R  =  ±  0.025 
or,  the  percentage  error  in  e  =  12.5%.  This  error  would  be  reduced  if 
a  higher  probe  temperature  (T  )  or  a  higher  critical  temperature  (T  ) 
were  used. 

An  analysis  of  the  error  in  computing  Rc and  t  is  more  difficult 
since  the  variables  in  these  quantities  change  with  time.  It  was 
estimated  that  at  short  times  the  error  in  both  Rc  and  x  could  be  of  the 
order  of  10%.  However,  this  error  decreases  significantly  as  time  is 
increased. 
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TABLE       II 


VALUES    CF   e 

LOCATION(R) 


FOD    NONDIMENSIONAL    TIME(X) 
tGNE-DIMENSIGNAL    MODEL. 


AND    NONDIMENSIOMAL 


R 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 
2.2 
2.4 
2.6 
2.8 

3.0 
3.2 
3.4 
3.6 
3.8 

4.0 

5.0 
6.0 
7.0 
8.0 
9.0 
10.0 


•^=0.1  —0.9 
R=1.0  —  10.0 
(3  =  0.0 


0.1 

1.000 
0.601 
0.316 
0.143 
0.055 


0.2 

1  .000 
0.691 
0.450 
0.274 
0.156 


0.3 


1. 

0. 
0. 
0. 
0. 


000 
733 
519 
352 
229 


0.4 

1.000 
0.758 
0.562 
0.405 
0.282 


0.5 

1.000 
0.776 
0.592 
0.443 
0.323 


0.6 

1.000 
0.789 
0.615 
0.472 
0.356 


0.7 

1.000 
0.800 
0.634 
0.496 
0.382 


0.8 

1.000 
0.808 
0.648 
0.515 
0.405 


0.9 

1.000 
0.815 
0.661 

0.532 
0.424 


0.018    0.082    0.142     0.190    0.230  0.263  0.290    0.313  0.333 

0.005    0.039    0.083    0.124    0.159  0.189  0.216    0.23^  0.259 

0.001    0.017    0.046    0.077    0.107  0.134  0.158    0.179  0.199 

0.000    0.007    0.024    0.047    0.070  0.092  0.113    0.132  0.150 

0.002    0.012    0.027    0.044  0.062  0.079    0.096  0.111 

0.001    0.006    0.015    0.027  0.040  0.054    0.068  0.081 

0.000    0.002     0.008    0.016  0.025  0.036    0.047  0.058 

0.001    0.004    0.009  0.016  0.024    0.032  0.041 

0.000    0.002    0.005  0.009  0.015    0.021  0.023 

0.001    0.002  0.005  0.009    0.014  0.019 

0.000    0.001  0.003  0.006    0.009  0.013 

0.000  0.000  0.000    0.001  0.001 

0.000  0.000 
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TA3LE   III 

VALUES  OF  G  FOR  NONOIMENSIONAL  TIME(r)  AND  NONOIMENSIONAL 
LOCATION(R)  f ONE- D I  MENS  I ONAL  MODEL. 

t»l.  0—9.0 
R=1.0  —  10.0 
§=0.0 


R 
1.2 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

0.821 

0.854 

0.8  69 

0.379 

0.885 

0.890 

0.894 

0.898 

0.900 

1.4 

0.672 

0.732 

0.  766 

0.777 

0.789 

0.798 

0.805 

0.811 

0.816 

1.6 

0.546 

0.627 

0.665 

0.689 

0.706 

0.718 

0.723 

0.736 

0.743 

1.8 

0.440 

0.536 

0.583 

0.612 

0.633 

0.648 

0.661 

0.671 

0.679 

2.0 

0.351 

0.458 

0.511 

0.544 

0.568 

0.586 

0.601 

0.612 

0.632 

2.2 

0.277 

0.389 

0.447 

0.434 

0.510 

0.531 

0.547 

0.560 

0.571 

2.4 

0.216 

0.330 

0.390 

0.430 

0.45  9 

0.481 

0.493 

0.513 

0.525 

2.6 

0.166 

0.278 

0.340 

0.332 

0.412 

0.435 

0.454 

0.470 

0.483 

2.8 

0.126 

0.233 

0.296 

0.338 

0.370 

0.394 

D.414 

0.430 

0.444 

3.0 

0.094 

0.1°4 

0.256 

0.299 

0.332 

0.357 

0.377 

0.394 

0.409 

3.2 

0.069 

0.161 

0.221 

0.264 

0.297 

0.32  2 

0.344 

0.361 

0.377 

3.4 

0.050 

0.132 

0.190 

0.233 

0.265 

0.292 

0.313 

0.331 

0.347 

3.6 

0.036 

0.108 

0.163 

0.205 

0.237 

0.263 

0.285 

0.303 

0.319 

3.8 

0.025 

0.083 

0.139 

0.179 

0.211 

0.237 

0.259 

0.277 

0.292 

4.0 

0.017 

0.071 

0.118 

0.157 

0.188 

0.213 

0.235 

0.253 

0.269 

5.0 

0.002 

0.021 

0.  049 

0.  076 

0.101 

0.122 

0.  142 

0.159 

0.174 

6.0 

0.000 

0.035 

0.018 

0.034 

0.0^1 

0.067 

0.  082 

0.096 

0.110 

7.0 

0.001 

0.005 

0.G13 

0.023 

0.034 

0.045 

0.056 

0.067 

8.0 

0.000 

0.001 

0.005 

0.010 

0.016 

0.024 

0.031 

0.039 

9.0 

0.000 

0.001 

0.004 

0.007 

0.  012 

0.016 

0.022 

10.0  0.000    0.001    0.003    0.005    0.008    0.012 
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TABLE  IV 


VALUES  OF©  FOR  NONDI  MENS  I ONAL  TIME(-c)  AMD 
NCNDIMENSIONAL  LOCATION(R),  ONE-DIMENS I ONAL 
MODEL. 


r.  =10.0  —loo, 

,0 

R=1.0 

— 10.0 

$=0.0 

T 
R 
1.0 

10.0 

20.0 

30.0 

40.0 

50.0 

70.0 

100.0 

1.000 

1.000 

1.000 

1.000 

1.00  0 

1.000 

1.000 

1.2 

0.903 

0.916 

0.922 

0.926 

0.929 

0.933 

0.937 

1.4 

0.820 

0.845 

0.857 

0.864 

0.869 

0.877 

0.8  84 

1.6 

0.749 

0.783 

0.800 

0.810 

0.818 

0.828 

0.838 

1.8 

0.687 

0.729 

0.750 

0.763 

0.772 

0.785 

0.797 

2.0 

0.631 

0.681 

0.705 

0.720 

0.731 

0.746 

0.760 

2.2 

0.581 

0  . 6?  7 

0.664 

0.68  2 

0.694 

0.711 

0.728 

2.4 

0.536 

0.  597 

0.628 

0.647 

0.660 

0.6"70 

0.698 

2.6 

0.494 

0.361 

0.594 

0.615 

0.629 

0.650 

0.670 

2.8 

0.457 

0.529 

0.563 

0.585 

0.601 

0.623 

0.644 

3.0 

0.422 

0.497 

0.534 

0.557 

0.574 

0.  598 

0.621 

3.2 

0.390 

0.46« 

0.507 

0.532 

0.550 

0.575 

0.598 

3.4 

0.360 

0.441 

0.432 

0.503 

0.5  26 

0.5^3 

0.578 

3.6 

0.333 

0.4ito 

0.458 

0.485 

0.505 

0.53  2 

0.558 

3.8 

0.307 

0.392 

0.436 

0.464 

0.484 

0.512 

0.539 

4.0 

0.283 

0.370 

0.415 

0.444 

0.465 

0.494 

0.522 

5.0 

0.188 

0.277 

0.325 

0.358 

0.381 

0.414 

0.446 

6.0 

0.122 

0.207 

0.256 

0.290 

0.315 

0.350 

0.385 

7.0 

0.077 

0.153 

0.201 

0.235 

0.260 

0.297 

0.334 

8.0 

0.047 

0.112 

0.157 

0.  1^0 

0.216 

0.253 

0.291 

9.0 

■   0.027 

0.081 

0.122 

0.153 

•0.178 

0.215 

0.254 

10.0 

0.015 

0.057 

0.094 

0.  123 

0.147 

0.183 

0.222 

20.0 

0.000 

0.001 

0.003 

0.009 

0.016 

0.031 

0.053 

30.0 

0.000 

0.000 

0.000 

0.001 

0.003 

0.010 

40.0 

0.000 

0.000 

0.002 

50.0 

0.000 
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VALUES    DF   e 

LOCATION(R) 


TABLE       V 

FOR    MOMDTMENSIONAL    TIMECC) 
, ONE-DIMENSIONAL    MODEL. 


AND    NONDIMENSIONAL 


"C  =0.1 

R=1.0 

(^  =  0.  1 

—  0.9 
-10.0 

X 
R 

1.0 
1.2 
1.4 
1.6 
1.8 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1. 
0. 
0. 
0. 
0. 

000 
,599 
,314 
,14  2 

055 

1.000 
0.688 
0.447 
0.272 
0.154 

1  .000 
0.729 
0.514 
0.348 
0.225 

1.000 
0.754 
0.555 
0.398 
0.277 

1.000 
0.770 
0.5  84 
0.435 
0.316 

1.000 
0.733 
0.606 
0.463 
0.347 

1.000 
0.793 
0.623 
0.485 
0.372 

1.000 
0.800 
0.637 
0.503 
0.392 

1.000 
0.807 
0.649 
0.518 
0.410 

2.0 
2.2 
2.4 
2.6 

2.8 

0, 
0. 
0. 
0. 

,018 

005 

,001 

,000 

0.081 
0.03<> 
0.017 
0.007 
0.003 

0.139 
0.092 
0.046 
0.024 
0.012 

0.186 
0.121 
0.  076 
0.045 
0.026 

0.224 

0.155 
0.104 

0.043 

0.255 
0.183 
0.129 
0.089 
0.059 

0.2  31 
0.203 
0.  151 
0.  103 
0.076 

0.302 
0.229 
0.171 
0.126 
0.091 

0.321 
0.248 
0.  189 
0.  142 
0.105 

3.0 
3.2 
3.4 
3.6 

3.8 

0.001 
0.000 

0.006 
0.002 
0.001 
0.000 

0.014 
0.008 
0.004 
0.002 
0.001 

0.026 

0.015 
0.  009 
0.005 
0.003 

0.039 
0.025 
0.015 
0.009 
0.005 

0.052 
0.035 
0.023 
0.014 
0.009 

0.064 
0.045 
0.031 
0.020 
0.013 

0.077 
0.055 
0.039 
0.027 
0.013 

4.0 
5.0 
6.0 
7.0 
8.0 
9.0 
10.0 

0.000 

0.001 
C.000 

0.003 
0.000 

0.005 
0.000 

0.009 
0.001 
0.003 

0.012 
0.001 
0.000 
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TABLE   VI 

VALUES  OF©  FOP  NOND IMEMSI ONAL  TIME  CO  AND  NONDIMENS IONAL 
LOCATION(R) , ONE-DIMENSIONAL  MODEL. 

T  =  1.0  —9.0 
R=l. 0—10.0 
£  =  0.1 

X         1-0         2.0         3.0         4.0         5.0         6.0  7.0         8.0         9.0 
R 

1.0       1.000    1.000    1.000    1.000    1.000    1.000  1.000    1.000    1.000 

1.2       0.812    0.841    0.853    0.860    0.864    0.867  0.869    0.«70    0.872 

1.4      0.658    0.710    0.732    0.744    0.752    0.757  0.761    0.764    0.766 

1.6       0.531     0.601    0.631    0.648    0.658    0.666  0.671    0.675    0.678 

1.8       0.425    0.508    0.545     0.565    0.579    0.588  0.594    0.599    0.603 

2.0       0.337    0.429    0.471    0.494    0.510    0.520  0.528    0.533    0.538 

2.2       0.265    0.362    0.407    0.433    0.450    0.462  0.470    0.476    0.481 

2.4       0.205    0.304    0.351    0.379    0.397    0.410  0.419    0.426    0.432 

2.6       0.157    0.254    0.303    0.332    0.351     0.365  0.375    0.382    0.388 

2.8       0.119    0.211    0.260    0.291    0.311    0.325  0.335    0.343    0.349 

3.0       0.089    0.175    0.224    0.254    0.275    0.239  0.300    C.308    0.314 

3.2       0.065    0.144    0.191     0.222    0.243    0.257  0.269    0.277    0.283 

3.4       0.047    0.118    0.163    0.193    0.214    0.229  0.240    0.256 

3.6       0.034    0.096    0.139    0.168    0.189    0.2C4  0.215    C.224    0.231 

3.8       0.024    0.077    0.118    0.146    0.168    0.181  0.193    0.201    0.208 

4.0       0.016    0.062    0.099    0.126    0.146    0.161  0.172    0.181     0.188 

5.0       0.002    0.019    0.040    0.059    0.075    0.038  0.093    0.106    0.112 

6.0       0.000    0.005    0.014    0.026    0.036    0.046  0.054    0.061    0.066 

7.0                        0.001    0.005    0.010    0.017    0.023  0.029    0.034    0.038 

8.0                        0.000    0.001    0.004    0.007    0.011  0.015    0.018    0.021 

9.0                                       0.000    0.001    0.003    0.005  0.007    0.009    0.012 

10.0                                                          0.000  .0.001    0.002  0.003    0.005    0.006 
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TABLE       VII 


VALUES    OFe    FOR    NOND I  MENS  I ONAL    TIME(-c)    AND 
NONDIMENSICNAL    LOCATION(R),     ONE-DIMENSIONAL 
MODEL. 

"C  =  10.( 
R  =  l  .0 
£  =  0.1 

)  -100. 
—  50.0 

,0 

R 

1.0 

1.2 

1.4 

1.6 

1.8 

10.0 

20.0 

30.0 

40.0 

50.0 

70.0 

100.0 

l.COO 
0.872 
0.768 
0.680 
0.605 

1  .000 
0.876 
0.774 
0.689 
0.616 

1.000 
0.876 
0.775 
0.690 
0.618 

1  .000 
0.876 
0.775 
0.690 
0.  618 

1.000 
0.876 
0.775 
0.691 
0.618 

1.000 
0.877 
0.775 
0.691 
0.618 

1.000 
0.877 
0.775 
0.691 
0.618 

2.0 
2.2 
2.4 
2.6 

2.8 

0.541 
0.485 
0.436 
0.392 
0.354 

0.553 

0.499 
0.451 
0.409 
0.371 

0.556 
0.501 
0.454 
0.412 
0.375 

0.556 
0.502 
0.454 
0.413 
0.375 

0.556 
0.502 
0.455 

0.413 
O.S"^ 

0.556 
0.502 
0.455 
0.413 
0.376 

0.556 
0.502 
0.455 
0.413 
0.376 

3.0 
3.2 
3.4 

3.6 

3.8 

0.31° 
0.288 
0.261 
0.236 
0.214 

0.338 
0.308 
0.281 
0.257 
0.235 

0.341 
0.311 
0.285 
0.261 
0.239 

0.342 
0.  312 
0.236 
0.261 
0.240 

0.042 

0.312 
0.286 
0.262 
0.240 

0.342 
0.313 
0.2  86 
0.262 
0.240 

0.342 
0.313 
0.286 
0.262 
0.240 

4.0 
5.0 
6.0 
7.0 
8.0 
9.0 

0.193 
0.118 
0.071 
0.042 
0.024 
0.014 

0.215 
0.  140 
0.092 
0.061 
0.040 
0.02£ 

0.219 
0.144 

0.0G6 
0.065 
0.044 
0.030 

0.220 
0.145 
0.09  7 
0.  066 
0.045 
0.031 

0.22.0 
0.145 
0.097 
0.066 
0.0^5 
0.031 

0.220 
0.145 
C.0°8 
0.066 
0.045 
0.031 

0.220 
0.145 
O.OQ*3 
0.066 
0.045 
0.031 

10.0 
20.0 
30.0 
40.0 
50.0 

0.C07 
0.000 

0.017 
0.000 

0.020 

0.000 

0.021 
0.001 
0.000 

0.022 
0.001 
0.000 

0.022 
0.001 
0.000 

0.022 
0.001 
0.000 
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TABLE       VIII 

VALUES    OF©    FOR    NONDI MENSI ONAL    TIME(r)     AND    NONDIMENS I ONAL 
LOCATION(R) , ONE-DIMENSIONAL    MODEL. 

"C=0.1  -0.9 
R=l. 0  —  10.0 

£  =  0.5 

T        O.L         0.2         0.3         0.4         0.5         0.6  0.7         0.8         0.9 
R 

1.0       1.000    1.000    1.000    1.000    l.OCO    1.000  1.000    1.000    1.000 

1.2       0.592    0.676    0.713    0.735    0.749    0.759  0.766    0.772    0.777 

1.4      0.308    0.433    0.494    0.530    0.555    0.572  0.585    0.596    0.604 

1.6       0.139    0.261    0.330    0.374    0.404    0.427  0.444    0.457    0.468 

1.8      0.053    0.147    0.211    0.256    0.289    0.314  0.333    0.348    0.361 

2.0      0.017    0.076    0.129    0.170    0.202    0.227  0.247    0.263    0.276 

2.2       0.005    0.037    0.075    0.110    0.138    0.161  0.180    0.196    0.20° 

2.4      0.001    0.016    0.042    0.068    0.091    0.112  0.129    0.144    0.157 

2.6       0.000    0.007    0.022    0.041    0.059    0.076  0.091    0.104    0.116 

2.8                      0.002    0.011    0.023    0.037    0.050  0.063    0.074    0.085 

3.0                      0.001    0.005    0.013    0.022    0.033  0.043    0.052    0.061 

3.2                      0.000    0.002    0.007    0.013    0.021  0.028    0.036    3.043 

3.4                                       0.001    0.003    0.007    0.013  0.018    0.024    0.030 

3.6                                       0.000    0.002    0.004    0.008  0.012    0.016    3.021 

3.8                                                        0.001    0.002    0.004  0.007    0.010    0.014 

4.0                                                      0.000    0.001    0.002  0.004    0.007    0.009 

5.0                                                                     0.000    0.000  0.000    0.001    0.001 

6.0  0.000    0.000 
7.0 
8.0 
9.0 
10.0 
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TABLE   IX 


VALUES  OF  e  FOR  NOND IM ENS IONAL  TIMECC)  AND  NOND I  MENS  I ONAL 
LOCATION (R) , ONE-DIME NS IONAL  MOOEL. 

"Csl.  0  —  9.0 

R=1.0  —  10.0 
£=0.5 


R 

1.0 
1.2 
1.4 

1.6 
1.8 

2.0 
2.2 
2.4 
2.6 

2.8 

3.0 

3.2 
3.4 
3.6 
3.8 

4.0 
5.0 
6.0 
7.0 
8.0 
9.0 
10.0 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0        7.0 


8.0 


9.0 


1.000    1.000  1.000    1.000  1.000    1.000    1.000    1.000  1.000 

0.780    0.797  0.901    0.803  0.804    0.804    C.804    0.834  0.804 

0.610    0.640  0.648    0.652  0.653    0.653    £.654    0.654  0.654 

0.477    0.517  0.528    0.532  0.534    0.535    B.536    0.536  0.536 

0.371    0.419  0.432    0.438  0.440    0.441    £.441    0.442  0.442 


0.287 
0.220 
0.167 
0.126 
0.094 

0.069 
0.050 
0.036 
0.025 
0.018 


0.340 
0.276 
0.223 
0.181 
0.146 


0.355 

0.2C3 
0.241 
0.199 
0.165 


0 
0 
0 
0 
0 


lis 

094 
075 
060 
047 


C 

0 
0 
0 
0 


136 
112 
0^2 
078 
063 


0.361 
0.29Q 
0.248 
0.207 
0.172 

0.144 

0.120 
0.100 
0.083 
0.070 


0 
0 
0 
0 
0 


364 

30  2 
251 
210 
176 


0.147 
0.123 
0.104 
0.087 
0.0^3 


0.365 
0.302 
0.253 
0.211 
0.177 

0.149 
0.125 
0.105 
0.089 
0.075 


0.366 
9.  304 
0.254- 
0.212 
I).  178 

D.  149 
D.126 
tt.  106 

D.090 
0.076 


0.366 
0.304 
0.254 
0.213 
0.176 

0.150 

0.126 
0.106 
0.0°0 
0.076 


0.366 
0.304 
0.2  54 
0.213 
0.179 

0.150 
0.126 
0.107 
0.090 
0.076 


0.017    0.037    0.051     0.058  0.061  0.063    0.064  0.064  0.065 

0.001     0.01C    0.018    0.023  0.026  0.027    0.028  0.028  0.029 

0.000    0.002    0.006    0.009  0.011  0.012    0.012  0.013  0.013 

0.000    0.002    0.003  0.004  0.005    0.005  0.006  0.006 

0.000    0.001  0.002  0.002   0.002  0.002  0.003 

0.000  0.001  0.001    0.001  0.001  0.001 

0.000  0.000    0.000  0.000  0.001 
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TABLE 


VALUES    OF   e    FOR 
LOCATION(R)  tONE- 


NCNDIMENSIONAL    TIME 
•DIMENSIONAL    MODEL. 

11=0.1  —0.9 

R=l.  0—10.0 
£=1.0 


(-c)     AND    NONDIMENSICNAL 


R 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 
2.2 
2.4 
2.6 
2.8 

3.0 
3.2 
3.4 
3.6 
3.8 

4.0 

5.0 
6.0 
7.0 
8.0 
9.0 
10.0 


0.1 

1.000 
0.583 
0.301 
0.135 
0.052 

0.017 

0.005 
0.001 
0.000 


0.2 

1.000 
0.662 
0.417 
0.248 
0.138 

0.071 
0.034 
0.015 
0.006 
0.002 

0.001 
0.000 


1 
0 
0 
0 
0 


0.3 

000 
695 
471 
309 
195 


0 

0 
0 


118 
069 
038 
0.020 
0.010 


0.005 
0.002 
0.001 
0.000 


0.4 

1.000 
0.713 
0.501 
0.346 
0.233 

0.153 
0.097 
0.  059 
0.035 
0.020 

0.011 
0.006 
0.003 
0.001 
0.001 

0.000 


0.5 

1.000 
0.724 
0.521 
0.370 
0.259 


0 
0 
0 

c 

0 


178 
120 
078 
050 


0.019 
0.011 
0.006 
0.003 
0.002 


0 
0 


001 
000 


0.6 

1.000 
0.72  2 
0.534 
0.388 
0.2  78 


0 
0 
0 

r\ 


0 

o, 

0 
0 
0 

0 
0 


197 
137 
094 

n  ^  o. 


031  0.041 


026 
01^ 
010 
006 
003 

002 
000 


0.7 

1.000 
0.737 
0.544 
0.400 
0.292 

0.211 

0.151 
0.106 
0.074 
0.050 

0.034 
0.022 
0.  014 
0.009 
0.005 


0, 
0, 


003 
000 


0.8 

1.000 
0.741 
0.551 
0.409 
0.303 


0. 
0. 
0. 
0. 
0. 

0. 
0. 
0. 
0. 
0. 

0. 

0. 


222 
162 

116 

083 
058 

040 
027 
018 
012 
003 

005 

000 


0.9 

1.000 
0.744 
0.55o 
0.416 
0.311 

0.231 
0.170 
0.125 
0.090 
0.065 

0.046 
0.032 
0.022 
0.015 
0.010 

0.007 
0.001 
0.000 
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TABLE   XI 


VALUES  CP  © 

LOCATION(R) 


FOR  NONDIMENSIONAL  TIME(X) 

,ONE-DIMEMSIONAL    MODEL. 


AND    NONDIMENSIONAL 


T  =  1.0 

R=1.0 

(2  =  1.0 

9.0 
10.0 

R 

1.0 

1.2 

1.4 

1.6 

1.8 

1.0 

2.0 

•: 

1.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

1.000 
0.746 
0.560 
0.422 
0.317 

1.000 
0.755 
0.575 
0.422 
0.341 

1. 
0. 
0. 
0. 
0. 

000 
756 

578 
445 
345 

1.000 
0.756 
0.579 
0.446 
0.346 

1.000 
0.756 
0.579 
0.446 
0.347 

1.000 
0.756 
0.579 
0.446 
0.347 

1.000 
0.757 
0.579 
0.446 
0.347 

1.000 
0.757 
0.579 
0.446 
0.347 

1.000 
0.75"? 
0.579 
0.446 
0.347 

2.0 
2.2 
2.4 
2.6 

2.8 

0.238 
0.177 
0.132 
0.097 
0.071 

0.264 
0.205 
0.160 
0.124 
0.097 

0. 

0. 
0. 
0. 
0. 

?69 

210 
166 
1  30 
102 

0.270 
0.212 
0.167 
0.131 
0.104 

0.270 
0.212 
0.167 

0.131 
0.  104 

0.270 
0.212 
0.167 
0.131 
0.104 

0.271 
0.212 
0.  167 
0.132 
0.  104 

0.27  1 
0.212 
0.167 
0.132 
0.104 

0.271 
0.212 
0.167 
0.132 
0.104 

3.0 
3.2 
3.4 
3.6 

3.8 

0.051 
0.036 
0.026 
0.018 
0.012 

0.075 

0.058 
0.045 
0.035 
0.027 

0. 
0. 
0. 
0. 
0. 

081 
064 
0  50 
040 
0  31 

0.082 
0.065 
0.052 
0.041 
0.033 

0.082 
0.065 
0.052 
0.041 
0.03  3 

0.082 
0.065 
0.052 
0.041 
0.033 

0.08  2 
0.  065 
0.052 
0.042 
0.033 

0.033 
0.066 
0.  052 
0.042 
0.033 

0.083 
0.066 
0.052 
0.042 
0.033 

4.0 
5.0 
6.0 
7.0 
8.0 
9.0 
10.0 

0.008 
0.001 
0.000 

0.021 
0.005 
0.001 
0.000 

0. 
0. 
0. 
0. 
0. 

025 
007 
0  02 
001 
000 

0.026 

0.008 
0.003 
0.001 
0.000 

0.026 
0.009 
0.003 
0.001 
0.000 

0.026 

0.009 
0.003 
0.001 
0.000 

0.026 
0.009 
0.003 
0.001 
0.000 

0.027 
0.009 
0.003 
0.001 
0.  000 

0.027 
0.009 
0.003 
0.001 
0.000 
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TABLE   XII 


VALUE?  OF  ©  FOR  NONDIMENS IONAL  TIME  CO  AND  NONDIMEN- 
SIONAL  LDCATION(R),  TWO-DIMENSIONAL  MODELITRUMP 
SOLUTION) 


"C=0.5 
=  1.0 
=0.0 

X=1P 


<§ 


50.0 
30.0 


X 

R 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 


0.5 

0.492 
0.216 
0.082 
0.027 
0.008 
0.002 
0.000 


1.0 

0.583 
0.330 
0.176 
087 
039 
016 
OOo 


0 
0 
0 
0 


0 

0. 

0 

0. 

0 

0 

0 


3.0 

696 
497 

355 
249 
170 
113 
074 


5.0       10.0       20.0       50.0      100.0 


0.741 
0.568 
0.440 
0.339 
0.258 
0.194 
0.142 


0.775 
0.622 
0.513 
0.420 
0.353 
0.286 
0.241 


0.791 
0.650 
0.550 
0.464 
0.403 
0.344 
0.305 


0.803 
0.671 
0.576 
0.499 
0.440 
0.390 
0.350 


0.830 
0.723 
0.633 
0.573 
0.512 
0.470 
0.425 


5.0 
6.0 

7.0 
8.0 


0.002    0.048  0.105  0.193    0.264  0.316    0.392 

0.000    0.017  0.050  0.123    0.202  0.258    0.327 

0.005  0.022  0.077    0.152  0.214    0.278 

0.002  0.010  0.048    0.117  0.179    0.238 


10.0 
14.0 
20.0 
30.0 


0.000    0.002    0.017    0.065    0.125    0.174 

0.000    0.002    0.015    0.055    0.0°? 

0.000    0.O01    0.012    0.035 

O.COO    0.001    0.006 
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TABLE   XIII 


VALUES  OF  e  FOR  MONO  I  MENS IONAL  TIME  CO  AND  NONDIMEN- 
SIONAL  LOCATION(R),  TWO-DIMENSIONAL  MCHEL(TRUMP 
SOLUTION) 


P.  =  1.0 
-c=0.5 
£=0.0 

X  =  5R 


30.0 
50.0 


X 
R 

1.5 
2.0 
2.5 
3.0 
3.5 
4.0 
4.5 


0.5 

0.491 
0.216 

0.083 
0.028 
0.008 
0.002 
0.000 


1.0 

0.582 
0.330 
0.176 
0.087 
0.040 
0.016 
0.008 


0 
0. 

0 
0 
0 

o. 

0 


3.0 

695 
496 
354 
247 
170 
112 
073 


5.0      10.0      20.0      50.0      100.0 


0.739 
0.  56  6 
0.433 
0.337 

0.256 
0.  192 
0.072 


0.770 
0.617 
0.506 
0.413 
0.343 
0.280 
0.141 


0.786 
0.643 
0.539 

0.453 
0.388 
0.332 
0.232 


0.796 
0.661 
0.563 
0.486 
0.424 
0.374 
0.290 


0.819 
0.699 
0.611 
0.541 
0.485 
0.437 
0.332 


5.0 
6.0 
7.0 
8.0 


0.002    0.047  0.103  0.188    0.251  0.297    0.361 

0.000    0.017  0.049  0.119    0.138  0.239    0.300 

0.005  0.021  0.073    0.141  0.196    0.253 

0.002  0.009  0.045    0.107  0.163    0.217 


10.0 
14.0 
20.0 
30.0 


0.000 


0.002    0.015 

0.05  7 

0.112 

0.158 

0.000    0.001 

0.013 

0.047 

0.034 

0.000 

0.001 

0.011 

0.032 

0.000 

0.001 

0.006 
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TABLE       XIV 

VALUES    OF  e    FOR    NONDIMENS IONAL    TIMECC)    AND    NONDIMEN- 
SIONAL    LOCATION(R),     TWO-DIMENSIONAL    MODEL(TRUMP 
SOLUTION) 


T=0.5  -50.0 
=  1.0  —  30.0 
=  0.0 

X=8.5R 


I 


X  0.5         1.0         3.0  5.0       10.0      20.0  50.0       100.0 
R 

1.5  0.485    0.569  0.664  0.700    0.724    0.736  0.745    0.772 

2.0  0.211    0.315  0.455  0.512    0.552    0.573  0.588    0.623 

2.5  0.080    0.165  0.313  0.330    0.431    0.456  0.477    0.531 

3.0  0.027    O.OPO  0.213  0.2G4    0.345    0.376  0.400    0.450 

3.5  0.008    0.036  0.142  0.211    0.274    0.310  0.339    0.403 

4.0  0.002    0.015  0.093  0.155    0.224    0.264  0.295    0.350 

4.5  0.000    0.006  0.060  0.112    0.179    0.223  0.258    0.322 

5.0  0.002  0.038  0.081    0.145    0.195  0.231    0.286 

6.0  0.000  0.013  0.037    0.088    0.143  0.184    0.241 

7.0  0.004  0.015    0.053    0.107  0.150    0.201 

8.0  0.001  0.006    0.032    0.031  0.124    0.174 


0.002 
0.000 

0.038 
0.013 
0.004 
0.001 

0.  081 
0.037 
0.015 
0.006 

0.145 
0.088 
0.053 
0.032 

0.195 
0.  143 
0.107 

0.031 

0.000 

0.001 

o.ooc 

0.010 
0.002 
0.000 

0.044 
0.017 
0.001 
0.000 

10.0  0.000    0.001    0.010  0.044    0.085    0.128 

14.0  O.OOC    0.002  0.017    0.042    0.073 

20.0  0.000  0.001    0.007    0.027 

30.0  0.000    0.000    0.004 
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TABLE       XV 


VALUES  OF  e  FOR  NOND I MENSI ONAL  TIME 
SIGNAL  LOCATION(R),  TWO-DIMENSIONAL 
SOLUTION) 

-c=0.5—  50.0 
P=l  .0  —  30.0 

e=o.o 

X=11.5R 


(x)     AND    NONDIMEN- 
MOOEHTRUMP 


X 
R 

1.5 
2.0 
2.5 
3.0 
4.0 
4.5 

5.0 
6.0 
7.0 
8.0 

10.0 
14.0 
20.0 
30.0 


0.5 


1.0 


3.0 


5.0       10.0      20.0      50.0       100.0 


0.010    0.029    0.094  0.137    0.183    0.223  0.241    0.293 

0.006    0.022    0.083  0.124    0.172    0.214  3.231    0.277 

0.003    0.014    0.066  0.106    0.153    0.193  0.214    0.266 

0.001    0.008    0.050  0.086    0.134    0.178  0.197    0.245 

0.002    0.024  0.051    0.095    0.142  0.165    0.212 

0.001    0.016  0.037    0.075    0.125  0.150    0.199 

0.000    0.010  0.027    0.064  0.111    0.138  0.184 

0.004  0.012    0.040  0.084    0.118  0.160 

0.001  0.005    0.024  0.065    0.098  0.140 

0.000  0.002    0.014  0.049    0.083  0.123 

0.000    0.004    0.027  0.058    0.095 

0.000    0.007  0.027    0.061 

0.001  0.004    0.022 

0.000  0.000    0.003 
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TABLE       XVI 


VALUES    OF©    FOR    NONDIMENSIONAL    TIME(-c)     AND    NONDIMEN- 
SIONAL    LGCATIONIR) i     TWO-DIMENSIONAL    MODtL(TRUMP 
SOLUTION) 


"C  =0.5  —  50.0 
R=1.0  ^30.0 
£  =  0.0 
X=13.5R 


X 

R 

1.5 

2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 
6.0 
7.0 
8.0 

10.0 
14.0 
20.0 
30.0 


0.5 


1.0 


3.0 


5.0       10.0       20.0      50.0       100.0 


0.000    0.001    0.013    0.030    0.060    0.102    0.113    0.166 
0.000 


0. 

,012 

0. 

,027 

0. 

,058 

0.101 

0. 

116 

0.156 

0. 

,010 

0, 

,024 

0, 

,053 

0.095 

0, 

112 

0.153 

0, 

,008 

0. 

,020 

0. 

,048 

0.090 

0. 

103 

0.147 

0. 

,006 

0. 

,016 

0. 

.042 

0.083 

0. 

102 

0.146 

0, 

,004 

0. 

,012 

0. 

,036 

0.077 

0. 

097 

0.136 

0, 

,003 

0, 

,009 

0. 

,0  30 

0.069 

0. 

092 

0.  133 

0, 

,002 

0. 

,006 

0, 

,025 

0.063 

0. 

087 

0.124 

0. 

.001 

0. 

,003 

0, 

,016 

0.050 

0. 

076 

0.113 

0. 

,000 

0. 

,001 

0. 

,010 

0.039 

0. 

067 

0.101 

0. 

,001 

0. 

,006 

0.031 

0. 

058 

0.091 

0. 

,000 

0. 

,002 

0.017 

0. 

042 

0.073 

0. 

,001 

0.004 

0. 

013 

0.041 

0. 

,000 

0.003 
0.000 

0. 
0. 

004 
000 

0.020 
0.0^3 
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TABLE  XVII.   EXPERIMENTAL  DATA 


RUN  NO. 

X 

■p(cm) 

T'cm) 

Xco 

\Co 

Tp(°o 

ec 

refp 

ref 

R 

i 

J.o 

onz 

OAQ 

22. c 

7A.<\ 

46.5" 

0.28Z. 

n.cz. 

in 

2.^3 

10-  ii-m 

a.o 

0-85 

O.S3 

22. G 

2.%* 

485" 

o.Z&Z 

T.fea 

4ii 

2.51 

C  v~a_el ;  c^l) 

4.  o 

1.  oc 

o.<£3 

22  C 

2.1.1 

4«-S" 

o-ZSZ. 

T.CZ 

411 

3o3 

X=iR 

G.o 

l.i6 

o.lE 

22.C 

Z1.1 

48. 5 

O.  Z<S2 

7.C2 

4n 

341 

8o 

1.24 

o.78 

a2.c 

21.1 

48.  S" 

0.3.S2. 

7.ca 

4.11 

3.1C 

JO.O 

1.32 

o.S3 

2  2.G 

2i-ci 

A«.S 

ozaa 

7  62. 

4n 

-3f<7 

z 

o.S" 

&.11 

0.45 

23.7 

w.i 

48.5" 

O.Z44 

T.3Z 

4.si 

2i4 

M-  S-lt 

1.5" 

o.fcl 

O.SI 

23.  <T 

2ll 

48. -^ 

o-  Z44 

132 

4-51 

2.4iT 

(\~cxdi<x.V) 

3.© 

o.l7 

o.Go 

23.<i 

2.1.1 

46. 5" 

6-244 

7.3Z1 

4-tii 

2.9/ 

fc^lfc 

4.o 

l.ol 
1./3 

o.G7 
oil 

23.9 

2S\  .7 

4S-5 

<>244 

1-32 

4-51 

3.Zo 

5.o 

23.  <? 

21.4 

-18.5" 

C-244 

1.3a 

-t.tfi 

3Ao 

n.o 

o.l^ 

Z3A 

23.7 

24. 1 

48.5" 

0-244 

132 

4.«51 

3.C3 

l-*2C 

o.TS 

21. 1 

^8.3" 

0-244 

1.32 

4-57 

3.18 

3 

2.o 

o.<?<* 

OG?- 

23.  o 

2.1.1 

4*.  <5" 

o.zm 

1.34 

4.SH 

241 

ll-  8-  7( 

4.o 

I./7 

O.T3 

23. o 

21.1 

48.  £" 

to. 2.1  l 

7-34 

4-51 

3.5o 

Cv-cva.oui ) 

G.o 

/.ze 

o.So 

23-© 

21.1 

48*5' 

O   Z.1I 

7.34 

4-51 

3.SS" 

^*i(^ 

S.o 

!.4o 

o.87 

23-o 

271 

48.£ 

o.zn 

1.34 

4S1 

4.2o 

(O.O 

1. SO 

0.7  3 

2.3.0 

21. 1 

48.5" 

o-21\ 

1.34 

4-51 

4.4^ 

A 

a.o 

0.<?<? 

O.C2, 

2.3.0 

21.1 

48.^ 

o.21\ 

T..3Z 

4-sn 

J.«?7 

it-<* -71 

4.o 

I./7 

o.l3 

JL3.o 

21.1 

46.^ 

o.2i\ 

1.32 

4-si 

J.  So 

^*~<xci  iul) 

G>.o 

(-32. 

O.S3 

23.  © 

21.1 

48.^ 

o-2lt 

7.32 

4si 

3.11 

*  =  ift 

g.o 

t-4-l 

o.88 

23.  o 

21.  1 

48.5* 

©-27V 

7.32 

4S1 

4.2.3 

(  o.o 

I.  So 

0.7  3 

23. o 

21.1 

48- 5 

OZT| 

1.32 

4- si 

4.4*i 

5" 

a.o 

o.7l 

O.  S7 

A4.o 

2.1.  1 

4ft  •£ 

0-24& 

7. U 

4si 

2.82. 

H-/1-  7  I 

4.o 

1.12. 

o.nz. 

A4.o 

21-1 

48.5 

o-24o 

7.11 

4.57 

3.4C 

(jrctU  O-O 

X^iR. 

91 


TABLE  XVII.  EXPERIMENTAL  DATA  (Cont'd) 


RUN  NO 

X 

£(cm) 

Hem) 

X<°c 

T/c) 

T/o 

ec 

(cm) 

ref 

(cm) 

R 

6 

2o 

o.^f 

o.S<i 

23C 

271 

48.5 

0.253 

T-Ll 

4.SH 

2. S3 

1 1  -  zo  -  7  t 

G.o 

1-28 

o.<&3 

23.  G 

21-1 

48.5 

02S3 

7.11 

4. si 

3  78 

Cv-oA^q^O 

S.o 

1-38 

O.Q1 

23.C 

21. 1 

485 

G-2.S3 

Til 

4. sn 

4.28 

*--  1R 

I  o.o 

MS 

0.15 

23.  G. 

21-  1 

48.5 

0-ZS3 

7.11 

4.5^ 

4.SG 

7 

2.o 

O.SG 

o.5S 

3.42 

21-1 

48.5 

O.Z35 

£.83 

A.S^ 

2.ni 

i\-2i  -71 

4-o 

I.09 

0-T3 

24.2 

21.  *? 

48.5 

0.235 

G.83 

4. si 

3.<r\ 

O^cl^O 

e.o 

1.3  7 

o.U 

14.2. 

2<<.1 

46.5 

O.  23b 

G-83 

4-  51 

4  42 

K--1R 

lo.o 

14  2. 

o.^S 

24.2 

211 

48.  s" 

o.  23S 

683 

4-51 

451 

£ 

\  .o 

O.&l 

o.4l 

2.4-2 

21.  9 

48.5 

0235 

&.fc3 

4sn 

1.17 

,  ii-  2.1-qi 

3.o 

o.f'i 

o.QG 

24.  2 

2<*1 

49-5 

0-Z3s 

6.-83 

4.5H 

3.2o 

(.V-<XcU<i-0 

5To 

in 

O.T8 

2.4.2 

2^1 

48.5 

o-235 

6.<&3 

4.51 

3.T7 

X=iR 

7.o 

V.32 

o.^e 

242 

211 

40-5 

0-235 

G.*3 

4  sn 

4  ISo 

f.o 

1-44 

o.  16 

24.2 

21.1 

43-  5 

o.23b 

6<33 

45T 

4-&2 

o.T4 

0.4^ 

9 

1.5 

242 

21- <? 

49.  5 

0-23S 

G.83 

4.57 

2  37 

/'-  2l-rl( 

3.S" 

V.o3 

0.6**1 

242 

23  <t 

48.5 

O.Z3  5 

G-83 

4  .si 

3. 3  I 

Oc*dU<x.f\ 

5.5 

I. IS 

o.  Tl 

24.2. 
^24. 2 

2.1.1 

25-1 

485 

0-235 

£.83 

4  ^'n 

3.8  2 

J^=1R 

1.^" 

1-31 

o.S8 

48-5 

O.  235 

<S&3 

4.51 

4-2£. 

lO 

i.o 
1.5 

0-G5 

O.A3 

25.2 

2.1. 'r 

48- S7 

o.Zoa 

6.83 

4. ST 

2.0S 

/I-2Z-71 

o.84 
1-14 

o.5l 

0.68 
O.T6 

2S2 

2<f  <i 

48.5 

o-Zol 

£-8  3 

4. ST 

2.4G 

^TCcli  fcM 

Z.o 

252 

23.1 

48.i»' 

O.2o2 

6. S3 

451 

2. To 

X=  \R 

3.o 
3.5 

25Z 
2  5  Z~ 

21  .*» 
211 

48. 5" 
48.5 

o.2o2 
o-2o2 

£83 
&.<33 

4.5H 

3.ZS 

4.51 

3.<2.S 

A-o 

l.o<t 

o.13 

2S2 

231 

49.5 

oz.oa 

CQ3 

4.si 

3.^1 

S".o 

12G 
1-32. 

o.©4. 
o.©8 

25  2. 

2S2 

21.1 

271 

48.5 

o."2.c2 

6.83 

451 

4.o5 

5".  5 

48.5" 

o-Z.o<L 

G<33 

-4. SI 

4  2fc 

Co 
7.o 

1.1A- 
1.3S 

o.83 
o13> 

252 

211 

48.5 

o  2o2 

68  3 

4. si 

4o  I 

25.  2 

21. 1 

48.5 

o.2o2 

683 

4.^q 

4.4G 
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TABLE  XVII.  EXPERIMENTAL  DATA  (Cont'd) 


RUN  NO. 

T 

Qcm) 

r>m) 

T0(°c) 

\(°C) 

Tp(°c) 

©c 

refD 

(cm) 

ref 

(cm) 

R 

io  (conO 

7.£T 

1.41 

o  77 

25\2 

27.3 

4S-^ 

0.202 

to. S3 

4.51 

4.15" 

8.0 

1.40 

o  .74 

25\2 

27.7 

4s.  5 

0.202 

G.83 

4. 31 

4.50 

9o 

i.  s2 

f.o2 

25  2 

27.1 

485 

0.202 

6. S3 

4. sn 

4.71 

9.3 

I. Go 

l.ol 

252 

2<m 

48.iT 

0.202 

6.S3 

451 

5.  \g 

/0.  O 

1.52 

l.©2- 

as.  2 

21 7 

485 

e>.2©2 

6.83 

451 

4  7  1 

il 

2.o 

6.71 

0.4G 

22. S 

27. 7 

48-5 

0  2§S 

3.4o 

2.2  \ 

2.23 

(i-  28-n 

4.o 

0.7/ 

0.57 

22.5 

211 

4S.S 

0.28-S 

3.4o 

2. 2.1 

2«5& 

(comcA) 

C  o 

l.ofc 

o.GG 

22.S- 

21.1 

49.5 

©-'2S-S 

3.4e> 

2.2  1 

i.  IS 

X-    1.R 

S.o 

l.oT 

o.C<j 

iZ.-S 

27-7 

43.  s- 

0.  28S- 

3.4ti 

2.21 

3.34 

1  O.O 

l.iE 

O.T3 

22. S" 

27-  7 

4S.S 

o-285 

3.4^ 

2.21 

3.5-0 

X  =   g.Sft 

Z.o 

o.&i 

O.4o 

a2.i> 

27-7 

7S.S" 

0.  z<as- 

3  4o 

2.2.1 

1-7/ 

4.0 

O&G* 

o.43 

22. t> 

2.7.7 

48-  5 

0.  28S 

3.40 

2.21 

2.o1 

C  o 

o.  g/ 

O.ST3 

22- iT 

2-7-1 

48-3" 

O.  28<^ 

3.40 

2.2.1 

2..S4 

9.o 

o.q\ 

O.S7 

22-5T 

Z7.7 

48$ 

o-Z8s 

3.40 

2.2.1 

2..S6 

J  o.o 

o.<n 

0  63 

22.  t> 

27.7 

48.-S 

o.ZSS' 

3.40 

2-21 

3.o2- 

y.  =    ior. 

2..o 

o.s-| 

0.3* 

Z2.S" 

2.7.7 

48- S 

o-28^ 

3.40 

Z.2.1 

1.57 

4.o 

e>.s"fe 

0.3& 

E^.tT 

21.7 

40.  s 

o-2BS 

3.40 

2.2.V 

\-1S 

£».  o 

o.t<s, 

o.43 

22.S" 

2.7.7 

40-S" 

o2SS 

3.4c 

2.2-* 

Z.ol 

S.o 

o.Tj 

0.4& 

22.5 

2.1.7 

48-  S 

o-ZSS" 

3.4o 

2.21 

2.23 

1  o.  o 

o.% 

0.4^ 

22-  £" 

21.7 

48- <S 

o-2S^ 

34o 

22.1 

2.37 

*-    W.SG. 

a.o 

o.33 

o.2o 

22. 5 

27.9 

4©.^ 

a  2.8  5 

3.4o 

22.1 

l.oo 

4.o 

o.33 

o-  2o 

22-5 

297 

4«^ 

a.2<8S" 

3.4o 

22.1 

l.Oo 

G.o 

o.33 

Cj.£o 

22.tr 

271 

48-3- 

o2Q5 

3.40 

Z.ZA 

lOo 

B.o 

o.  3  3 

0.  2o 

22.5 

27. ^ 

40-S 

o28S 

3.  4  0 

2.2.1 

l.Op 

(O.o 

o.  4( 

o.?L& 

22  .s 

21-7 

4&<S 

OffiS 

3.4o 

22. 1 

121 

■ 

. 
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a: 

a: 

o 

C 

UJ 

UJ 

I/O 

h- 

K 

_l 

5: 

• 

—z 
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cc 

T-* 

t-4»— 1 
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> 

UJ 

CO 

OC 
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5: 

wU_ 

LL 

UJ 

< 
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O 

1- 

2 

a: 

nc 

-z. 

O 

o 

UOh 

K 

>— < 

_l 

o 

Z>"~ 

>—i 

UL 

OH 

rxs 

:>: 

u_ 

a. 

u.      ►-. 

»— 4 

o 

o 

z-j 

-J 

o 

a: 

-JO 

a: 

c 

UJ 

•<»— |Q£ 

or 

UJ 

-J 

H 

Z  CO  UJ  r-lLU  O  CO 

CO 

3 

arzrisoo- 

•  s:< 

D  It 

0. 

LUUJOOQ-C^i 

^< 

s: 

t-S-J 

•  !Dr-4  2T< 

Mi- 

o 

Xm  II 

II    II 

ii  ii 

ll  UJ 

o 

ujQ<<aj  co^rzrco 

< 


o 


a 

LL 


<  CO 

UJ 

<  KH 

CC 

cc  < 

o  > 

u 

UJ 
CO  5! 

UJ  M 

3                                      f^r-4     rr,  in  o  (Nj  (nj  >0  mr^JOirir  j  r-ir^rvjr-rriirN  ao  •-<  H- 
_i  -o  0s  r-  o  c  r-  o  <o  p-  <^  <f  o  o  ^o  r-  co  •-<  co  m  tr  *-<  ro  (\j  ro  oo  r-  ro  o  .-i  r~-  vo 
<f^(^vt^VHCOOo-^^>^co^^oO'HOC>i^0^^c^ro.rococ>^LnCTsco 

>^^c^rO(Njrv|r^vOc\jc\Jo0^^^vtiX>roONii".(^cv'»^o^^(NjroirtCOf-4-j-oo  • 
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